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Apparatus for the Determination of Minor Components 
of a Gas Mixture 


Martin Shepherdt 


The apparatus is primarily for the determination of the components occurring in low 
partial pressures in a gas mixture, and which ordinarily cannot be determined with satis- 
factory accuracy because of the significant sorption and desorption of major components 


in the reagents intended for the removal of the minor ones. 


Small amounts of solid and 


liquid reagents are used in the apparatus, and combustions, as well as absoprptions, can be 


performed. 
for volumetric chemical gas analysis. 


1. Introduction 


There comes a time in the life of every exasperated 
vas analyst when he tires of cleaning liquid reagents 
from parts of the apparatus not intended to contain 
them, and the tendency of solutions to wander out 
of bounds prompts the desire to substitute 
nomadic solid reagents. And there are times when 
the analyst is discouraged by conventional analyses 
that have indicated minor constituents that did not 
exist in the samples analyzed, or lost small but 
significant portions of major constituents. Errors 
caused by the solubility of components in a reagent 
not intended for their removal have long been a 
problem in gas analysis. During 1938 the apparatus 
described here was designed to eliminate such errors 
of solubility, at least within the significance of the 
better volumetric systems. Prior to this, the idea 
had been used in determining carbon dioxide in air 
captured in the stratosphere [1],’ and for the deter- 
mination of ethylene oxide in carbon dioxide [2]. 
Part of this apparatus was illustrated and _ briefly 
mentioned in a paper [3] describing modifications of 
what was then our volumetric gas analysis apparatus. 

The determination of constituents of small partial 
pressures has been successfully done by Bone and 
Wheeler, whose apparatus seems to be one of the 
standard types used in England [4]. This appara- 
tus requires the addition of a small amount of fresh 
reagent for each absorption conducted, and while 
this minimizes the errors of solubility associated 
with the ordinary volumetric apparatus (usually 
referred to as an Orsat), this gain in accuracy is 
obtained by a considerable inconvenience in opera- 
tion. Analysts in general have shown too great 
an impatience ever to accept the Bone and Wheeler 
procedure. So they have put up with errors of 
solubility. 

Errors caused by solubility have long been recog- 
nized and even talked about, but it took a demon- 
stration like the cooperative analysis of a sample of 
natural gas by 30 laboratories working with Sub- 
committee D-3—VII (Analysis of Gaseous Fuels) of 
the American Society for Testing Materials to bring 
a full realization of the magnitude of these errors [5]. 


less 


* Deceased, September 17, 1953, 
Figures in brackets indicate the literature references at the end of this paper. 


| of oxygen in the sample. 





The apparatus can be made a separate unit or part of a standard apparatus 


The determination of oxygen was _ particularly 
interesting because there was no significant amount 
One- or two-tenths of a 
percent of oxygen is usually reported in the analysis 
of a natural gas, and the extent to which oxygen 
really does occur has long been a question. Of 294 
analyses reported, 124 (only 42 percent) correctly 
indicated no oxygen; 58 (about 20 percent) indicated 
0.1 percent of oxygen, which may be considered the 
average error of solubility that experience has shown 
may be expected; and 69 (about 23 percent) indi- 
cated 0.2 percent of oxygen, which should have been 
the maximum error caused by solubility alone 
However, 24 and 19 analyses indicated 0.3 and 0.4 
percent, respectively, and these values round off a 
probability curve. 

Errors of solubility can be minimized by using 
small amounts of solid reagents and liquid reagents 
in reaction tubes designed to confine them. The re- 
action tubes can be connected to the distributor of 
the gas analysis apparatus in such a way that the 
manipulations involved are more easily performed 
than in the case with the ordinary volumetric 
apparatus employing liquids in pipet. Actually, the 
reaction tubes form an integral part of the distribu- 
tor, and a reaction is conducted by simply passing 
gas from the buret through the tube and into a reser- 
voir over mercury, just as though the gas were being 
passed through the distributor itself. 

The apparatus can be made a separate unit or a 
part of the conventional gas analysis apparatus. 


2. The Separate Unit 


The separate unit is independent of the regular 
one ordinarily used for volumetric analysis. It is 
possible to do a complete absorption-combustion 
analysis with this unit, even with gases of average 
composition, but it is primarily intended for the 
determination of components of a mixture that occur 
in small amounts. If the type of analysis performed 
requires the determination of components present 
in both large and small amounts, the apparatus 
described in section 3 is the one to use. 

The various reaction tubes of this unit are per- 
manently sealed to the distributor, which may be 
made in 1, 2, or 3 sections, as desired. A natural 














SIDE VIEW 


Figure 1. 


division would be to include absorption tubes in one 
section (see fig. 1, stopeocks 1 to 5, inclusive), com- 
bustion tubes in a second (cocks 6 and 7), and the 
reservoir and connection for the pipet containing 
alkaline pyrogallo! solution in the third (cocks 8 
and 9). In arranging the section for absorptions, 
the particular needs of the work to be done will 
often indicate a simplification of the apparatus shown. 

If the distributor is made in a single section, the 
total length (60 em) makes it advisable to use 10- 
mm outside diameter capillary tubing for the hori- 
zontal member, and the vertical members are 8-mm 
outside diameter, and all bores are 2 mm. When 
the distributor is divided, the sections are connected 
with 12/2 interchangeable spherical joints, such as 
terminate the distributor at both ends. 

The buret used with this unit is not shown, but 
is just like the one previously described !6], complete 
with the special compensator-manometer, except 
that the top 45 ml of the buret itself is a pear-shaped 
bulb. The stem of the buret is graduated in the 
standard manner [6], but as the over-all length is 
the same and the volume only half that of the reg- 
ular buret, the stem may be marked in intervals of 
0.1 ml rather than 0.2 ml, which is advantageous 
for analysis of this kind, and the standard water 
jacket and frame may be used. 

As previously noted, the reaction tubes are in- 
tegral parts of the distributor and are designed to 
use small amounts of either solid or liquid reagents. 
They are connected by stopeocks with ol bored so 
that in one position the reaction tube is bypassed, 
and in another position gas from the buret enters 
the bottom of the reaction tube and leaves the top 
traveling left into the distributor on its way to the 
mercury-filled reservoir. The various tubes will be 
briefly described. 
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Separate unit. 


Reaction tubes A and B (fig. 1) are alike, and are 
for solid reagents. The ends of these tubes are 
closed with flat caps ground to fit the ground flat 
flanges at the top and bottom of the tubes. These 
joints may be sealed with a low-melting thermoplss- 
tic cement or with a heavy lubricant. The caps 
have glass spurs to hold retaining springs or rubber 
bands in place. These tubes are ordinarily used 
with a desiccant, such as phosphorous pentoxide on 
glass wool or calcium sulfate, in the event a dry 
analysis is done; or with Ascarite for the removal of 
carbon dioxide and other acid gases. 

Reaction tube C may be used with solid reagents, 
or the central tube may be used with a liquid reagent 
on glass wool while the outer tubes contain solid 
reagents. For example, a few drops of fuming 
sulfuric acid in the middle tube, and a packing of 
glass wool and Ascarite in the outer tubes, may be 
used to determine small amounts of unsaturated 
hvdrocarbons in fuel gases; or calcium sulfate in the 
outer tubes and Hopealite in the middle tube will 
serve to remove small amounts of carbon monoxide 
from the gas sample. When a liquid reagent is used 
in the central tube it is introduced through holes in 
the top cap and top male grinding of the tube, which 
may be made to coincide by rotating the cap. For 
convenience in filling, the open position is at the 
front of the tube, and the holes are drilled at a down- 
ward angle (fig. 1) to permit insertion of the nozzle 
of a medicine dropper, from which the reagent is 
transferred. The excess and spent reagents are 
removed through a similar arrangement at the bottom 
of the tube, again using the medicine dropper 
Dropping bottles may be used for the convenient 
storage of reagents. 

Tube D is for liquid reagents, and is like the middle 
tube of C except for internal diameter. Glass wool 
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lass tape is ordinarily used as a packing to hold 

Jutions. The diameter of this tube will permit 

other packing—a central glass rod wound with a 
elix of glass wicking, which snugly fills the annular 

ace, thus providing a relatively long path for the 
as to travel. 

Tube E is designed to permit coating the packing 
with fresh reagent several times during an absorption 
One charge of reagent may also last for quite a few 
wnalvses. The solution is placed in the lower 

servoir and forced up over the packing by com- 

essing a rubber bulb attached to the top outlet of 
this reservoir. The in-seal tube just above the 

servoir allows the solution to ascend the reaction 

be without entering the left-hand capillary to 
topeock 5 (which is closed during this procedure). 
The solution is brought to an etched mark on the 
neck just above the reservoir (before each gas volume 
s measured) by a screw clamp on the rubber buib 
attached at the top of the reservoir 

Tube F is for combustion with a catalyst heated by 
an outer tube furnace. Copper oxide is ordinarily 


uw 


ised 

Tube G is a combustion tube containing a platinum 
helix fastened to two platinum leads of much larger 
diameter than the helix Contact to power leads is 
made through two mercury-filled wells (W,, W.) in 
the top cap. This cap is fastened to the combustion 
tube with a high-melting thermoplastic cement, and 
is cooled with strong air jets during a combustion 
\ packing of fine quartz sand acts as a safe flashback 
trap. The combustible premixed with 
oxvgen before burning. 

Tube H is a mereury -filled reservoir, which con- 
fines the gases passed through the reaction tubes, 
and delivers them back to the buret For conven- 
ience in displacing gas from the reservoir, a ground 
spherical ball N fits against ground 45-degree seat 
\l. The seat is ground snugly into the bottom lead 
of the T-cock 8, presenting essentially no dead space 
The ball cuts off the flow of mercury as the gas is 
delivered from the reservoir with a volumetric repro- 
ducibility better than could be obtained by adjust- 
ment to an etched mark, and with far greater ease 

\ pipet of the distributor-tip type [7] is connected 
This furnishes the nitrogen necessary 
as a flushing, or sweeping, gas used to transfer the 
sample as required throughout the apparatus. 


Yases 


are 


to stopcock 9. 


3. Procedure for an Analysis With the 
Separate Unit 


3.1. Preparation of Transfer Gas 


The sweeping or transfer gas that may be needed 
s prepared by absorbing oxygen from air, and storing 
over the alkaline pyrogallol solution in the pipet 
connected to stopcock 9. The recently developed 
pyro soiution should be used [8]. The distributor, 
all reaction tubes, and the arm of the manometer 
connected to the distributor are adjusted to the 
relerence compensator pressure. 


3.2. Procedure for an Absorption With Solid Reagents 


The gas sample is measured and then transferred 
at the approximate rate of 100 ml/min through the 
first reaction tube (with all other tubes bypassed) 
to the reservoir H. The gas contained in the arm 
of the manometer connected to the distributor is 


displaced into the buret, and then transferred 
through the reaction tube. If it is important to 


keep the reagent at the exit of the reaction tube 
fresh, the residual yas is returned from the reservoir 
to the buret with the reaction tube bypassed, but if 
the condition of the reagent at the exit of the reaction 
is not important, the residual gas is returned through 
the tube The cvele Is repeated until all significant 
reaction has occurred as shown by no further con- 
traction. The pressure in the buret, d str butor, 
and reaction tube are adjusted to the reference 
pressure, and the volume of the residual 
measured 


yas 1s 


3.3. Procedure for Absorption With Liquid Reagents 


When the nitrogen is prepared at the beginning 
of the analysis, and all reaction tubes are filled with 
it, the liquid reagents to be used are introduced at 
the top of each of the reaction tubes designed to 
contain them. A very slow flow of nitrogen from 
the buret will prevent aiv from entering the reaction 
tube during introduction of the solution, but if there 
is any reason air contamination, a few 
final passages into the pyrogallol may be made 
The final nitrogen balance is made after the intro- 
duction of the liquid reagents. The actual procedure 
for absorptions is the same as that outlined for 
absorptions with solid reagents 


to suspect 


3.4. Procedure for Combustion 


The procedure for combustion with copper oxide 
does not differ from the routine ordinarily described 
Inasmuch as carbon dioxide, which may be produced, 
is strongly absorbed on a cooling catalyst, it is best 
to start and stop with the catalyst hot. The initial 
nitrogen balance should be made with this in mind 

The combustion over the hot platinum helix is apt 
to be more complete than the conventional slow 
combustion. The residual sample from the absorp- 
tion analysis is stored in the reservoir while the 
oxygen air) for the combustion is measured, 
These two gases are mixed by passing back and forth 
between reservoir and buret, and then the mixture 
is passed over the glowing platinum helix and thence 
through the flashback trap to the reservoir at the 
approximate rate of 20 ml/min. The gas in the arm 
of the manometer connected to the distributor is 
displaced and passed through the combustion tube 
After this, the residual mixture in the reservoir may 
be below the explosive limit, and this gas may be 
passed back over the platinum to the buret How- 
ever, the safety of this procedure should be estab- 
lished by experience. The cycle is continued until 
no further significant reaction occurs. Carbon diox- 


(or 








ide in the products of combustion is absorbed in 
Ascarite contained in one of the reaction tubes. In 
some analyses it is advantageous to absorb the carbon 
dioxide during the combustion procedures. For ex- 
ample, methane and ethane (and other hydrocarbons) 
caleulated from contraction and carbon dioxide pro- 
duced are seriously in error, whereas calculation from 
contraction, carbon dioxide produced, and oxygen 
consumed, will give more accurate results [9]. In 
such cases, the removal of carbon dioxide during the 
course of the combustion will result in a more rapid 
and possibly more complete combustion. 


3.5. General Procedure 


When the gas sample, or a part of it, has entered 
a reaction tube, some of the residue remaining in the 
tube is reactive with respect to the rest of the 
analysis. It is obvious that with each successive 


reaction the gas in formerly used tubes must enter | 


the tube used for each new reaction. Bringing all 
the gas into the reactive zone is quite easily done 
with this apparatus. In the more conventional 
volumetric apparatus, passing the gas into a suc- 
cession of absorption pipets or correcting for reactive 
residues by the calibrated distributor procedure [6] 
is considerably more trouble. 


The apparatus must be operated in an area of | 


reasonably constant temperature. One single reac- 
tion tube when filled with reagent contains about 
3.5 ml of gas, so that an error of about 0.01 ml/deg C 
change in temperature could be made in comparing 
volumes before and after an absorption. The error 
with a double or triple reaction tube would be almost 
proportionately greater. 
a complete analysis, using every reaction pipet 
shown, would be about 0.2 mli/deg C. This is too 
great, and temperature corrections should be applied 
if the laboratory is not suitably conditioned. The 
volumes of the various reaction tubes are easily 
determined by evacuation and subsequent filling 
from the buret. For more accurate work involving 


very small changes in volume, reaction tubes of | 


smaller volume than those shown should be used. 
As this apparatus is likely to be put to special uses, 
it should be tailored to fit. 


4. Reaction Tubes Adapted to Standard 
Volumetric Apparatus 


The average analytical demand may not require a 
complete separate unit of the kind just described. 
Instead, the standard volumetric apparatus will 
usually be employed for the analysis of most gas 
mixtures, but as most mixtures are apt to contain 
one or more components in small percentages, the 
reaction tubes might well be employed for part of 
the analysis. Fortunately, they are neatly adapt- 
able to the standard apparatus. This may be done 
with a set of special connectors and interchangeable 
spherical ground-glass joints. The reaction tubes 


may be connected across 2 T-cocks of a distributor, 
or better yet, to one cock of the type designed to pass 


The accumulative error for | 





gas into and out of a closed tube, or bypass the tub: 


Such units are commercially available. These reac- 
tion tubes, their connectors, and the distributor 
section designed for them are shown in figures 2, 3 
and 4. The connectors and the basic reaction tube 
for solids or liquids upon glass wool or tape are show: 
in figure 2. he tube itself is shown at A. With 
connector D it forms a single reaction tube that fits 
the outlets of the special distributor stopcock shown 
at G, figure 4. With connector C, two of A make a 
double reaction tube, and with connector B= and 
another of A added to the double tube, a triple 
reaction tube is assembled. The special tube for 
absorptions requiring several applications of fresh 
reagent is shown at F, figure 3; its connector is F of 
figure 2. The slow-combustion tube is shown at G, 
figure 3; its connector is G of figure 2. The copper- 
oxide tube is shown at H, figure 3; it needs no con- 
nector unless keyed into two T-cocks of a standard 
distributor, in which case two of connector H (fig. 2) 
will serve to bring it to the proper upward position 
The mercury-filled reservoir, K of figure 3 may be 
connected to a T-cock of the standard distributor if 
space permits, but the combustion pipet of the 
standard apparatus will serve for this reservoir, 
albeit not so neatly. Figure 4 shows a five-cock 
standard distributor with a two-cock section for 
reaction tubes. The special section is placed next 
to the buret (which is to the right but not shown). 
This arrangement serves for many types of analysis 
involving absorption and combustion. The five 
T-cocks will accommodate the slow combustion pipet 
and the two pairs of absorption pipets prescribed in 
the American Society for Testing Materials standard 
method for the analysis of natural gases and other 
fuel gases of similar type [10]. A copper-oxide tube 
may be connected at stopcock H. Stopecock G will 
accommodate any of the reaction tubes, including 
the triple reaction tube mentioned in the above 
ASTM standard for the determination of the small 
amounts of unsaturated hydrocarbons that may be 
in some natural gases or similar fuels. 

It is apparent that this type of apparatus can 
conveniently meet the needs of both special and 
regular problems, and perhaps occasionally stimulate 
the imagination of the more enthusiastic analysts. 
For example, the possibility of conducting an entire 
analysis with dry gases has been mentioned. If this 
is done, it will probably be best to dry the sample 
en route to the buret, although water vapor in the 
sample may be determined volumetrically if the 
walls of the buret are carefully equilibrated with 
respect to the sample before the initial measurement 
of the sample. Working without water has distinct 
advantages. The buret remains clean, and the posi- 
tion of the mercury meniscus is determined with less 
eye strain and perhaps greater precision than is the 
case when films or droplets of water are on the buret 
walls. There is no error caused by a change in the 
amount of water above the mercury in the buret or 
by solution of any gas in the water here or in the 
combustion pipet. (Water formed in the combustion 
pipet evaporates into the dried gas passed back into 
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Figure 3. Special reaction tubes and reservoir. 
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it during the combustion and is removed by the 
desiccant It should be noted that prolonged use 
of Asearite during dry analyses may remove so much 
water from it that the rate of absorption of carbon 
dioxide on this reagent is considerably reduced 
Although the likelihood of a general conversion to 
dry analyses is not great, every analyst should have 
the satisfving experience of performing a few. On 
the other hand, a wet analysis can be made much 
more dependably by interposing a saturator in the 
reaction tube nearest the buret, rather than the 
desiccant proposed for a dry analysis. One of the 
greatest and most frequently occurring errors in 
volumetric analysis is caused by failure to saturate 
properly each measured gas with water vapor 
Samples of compressed gases, and residues returned 
from desiccating reagents, are often not completely 
saturated. The time required for saturation may 
be verv long if water vapor must be supplied from 
liquid in the lower parts of the buret because the 
vapor travels slowly by diffusion and convection in | 
tubes of this size 


~ 


|4] George Lunge and H. R. 


{9} 
(10) 


Joseph R 
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Mismatch Errors in the Measurement of 
Ultrahigh-Frequency and Microwave Variable Attenuators 


R. W. Beatty 


Che mismatch error in the measurement of ultrahigh-frequency and microwave variable 


attenuators is analyzed and expressions are derived. 


mismatch error in measuring the difference 


An example is given to show that the 


in attenuation between two attenuators is less 


than the sum of the mismatch errors obtained when measuring each attenuator individualls 


1. Introduction 


The error resulting from mismatched generator 
nd detector sections in the measurement of the 
attenuation of a single attenuator is well known.! 
The corresponding error in the measurement of 
changes in attenuation is important in the calibration 
of variable attenuators and in the calibration of 
large attenuators, using known pads as gage blocks 
or fixed attenuation standards. 

The calibration of a variable attenuator consists 
in measuring the change in the insertion loss as the 
attenuator dial moves from a zero or reference 
position to another position that is marked or can 
be read on a scale. The change in the insertion loss 
equals the change in attenuation if the attenuator is 
placed in a refleetionless, or matched, system. 
There is always a degree of uncertainty regarding 
the match, depending upon the accuracy of the 
instruments used to indicate or recognize matched 
conditions and upon the reflections from connectors. 
For this reason the change in the insertion loss 
cannot be considered to be exactly equal to the 
change in attenuation, and the difference is called 
the mismatch error. 


2. Expression for Mismatch Error 


A change in a variable attenuator from the 
reference position to another position is equivalent 
to removing one attenuator and inserting another 
attenuator in the cireuit. The insertion in 
decibels, of the attenuator’? corresponding to the 
reference position of the variable attenuator is 


loss, 


(l- 


20 logio 


Silo) 
Sio(1 


SV) 
rl) 


L 


The voltage-reflection coefficieats Tg and Ty, refer 
to the generator and load, respectively, and are 
measured at the terminals where the attenuator is 
nserted. The scattering coefficients* S),, Sj, and 


C.G. Montgomery, Technique of microwave measurements, p. 824 (McGraw- 
Book Co., Inc., New York, N. Y., 1947 
see Appendix 
he seattering coefficients Sy, Sw, and S» are in effect defined by eq 12 and 
2. See also, C. G. Montgomery, R. H. Dicke, and E. M. Purcell, Prin 
of microwave circuits, p. 146-151 (McGraw-Hill Book Co., Inc., New 
wk, N. Y., 1948 


S» refer to the attenuator, corresponding to the 
reference or zero position of the variable attenuator 
The corresponding insertion loss for a different setting 
of the variable attenuator is 


( §:.1.)(1 —S.I..)—(S..)T 
20 logy, , al\ 22! = 2) Tals * (2) 
: S}.(1—Tl 1) 


I Bg 
The change in the insertion loss is 


; l 
AL=L’—L=20 logy Re 20 logo 


(1—S,,l,)(1 —S2lr.)—(S}.)'-l 
+20 logis —s as ==!) (3) 
(1—Sy Po) 1 —Sal)— SPT CP: 


AL= A’—A-+te, 
where rg A, and € in eq (4) correspond to the three 
terms in eq (5). The error ¢€ must be subtracted 
from the change in the insertion loss to obtain the 
change in attenuation, A’—A. 


The error term can also be written * 


iV) 


XY S31.) 
r.) (1 


( r 
») 

€ 20) log 0 - ’ : , 

a r; SoI",) 


where IT, and T, are the input-voltage reflection 
coefficients of the attenuator terminated in a load 
having a voltage reflection coefficient Ty. 

The mismatch error in the measurement of a 
single attenuator can be obtained as a special case 
of eq (5). The reference attenuator vanishes in 
this case, changing I, to T',, Sy, to unity, and S,, and 
S» to zero. Substituting these values for T, and 
Sx» into eq (5) yields 


(1—Ir°,l,)(1—S.I 1) 
€ 20 login — = ? 


(6 
F< ») 


which corresponds to eq (24), page 826 of the refer- 
ence of footnote 1. 


‘Substitute T= S)-+ 








3. Evaluation of the Mismatch Errors 


It is possible in principle to evaluate the mismatch 
error by measuring the voltage reflection coefficients 
I. Te. Ti. Ti. Sx, and Sj, and substituting them 
into eq (5). 

In many cases, the magnitudes of the reflection 
coefficients can be determined, but their phases 
cannot be conveniently determined because of the 
limitations of a particular messuring apparatus. 
Equation (5) can then be used to find the limits of 
the mismatch error, permitting the phases of the 
reflection coefficients to have all possible values. 
The limit of error can be expressed in the form 


(1 T ile 1 ~ + S21 1)) (7) 
(IF|PiTo)QF|Satil) * 


€1imie = 20 logio 


and the corresponding limit of error for single 
attenuators is 


(1 + Inre)a = Sx!) 


8 
1¥/lols! (8) 


€ limit 20 logio 


For exemple, in order to reduce mismatch errors in 
the calibration of attenuators, the magnitudes of 
I, and T,, are made as small as possible, and their 
probable amplitude is estimated from the accuracy 
of the apparatus used to recognize matched condi- 
tions and from the known connector characteristics. 
It is difficult to accurately determine the phases of 
these small reflection coefficients, and the mismatch 
error can generally be determined not exactly, but 
within limits. 

An example will illustrate the determination of 
mismatch error. If the voltage standing-wave 
ratios og, OL, 1, 1. Fx, and oy [e=(1+\T})/(1—|P})] 
corresponding to Tq, Ty), [Ty), \Pi!, |Sx, and |S: 
are 1.1, 1.1, 1.2, 1.5, 1.2, and 1.5, the limits of 
error for the initial attenuator calculated from eq (8) 
are approximately +0.095 decibel. The corre- 
sponding limits of error for the final attenuator are 
approximately +0.185 decibel. The limits of error 
for the change in attenuation calculated from eq (7) 
are approximately + 0.242 decibel. It is seen that 
the mismatch error for the change in attenuation is 
less than the sum of the mismatch errors in measuring 
each attenuator individually. 

It is often possible te determine the phases of 
r;, Ti, Sw, and S,, enabling one to calculate the 
mismatch error within narrower limits. The error 
in measuring voltage standing-wave ratios and phase 
angles produces an error in determining the limits of 
mismatch error but is in general a second-order 
effect that can be neglected. 





4. Appendix. Derivation of Expression fo: 
Insertion Loss 


The insertion loss of a four-terminal network j 
defined by the equation 


P, 

L=10 logy P? 
where P;, denotes the power delivered directly to a 
load, and P, denotes the power delivered to that 
load when the network is inserted between the load 
and the generator. It is assumed that the generator 
is stable and unaffected by changes in loading. 

The insertion loss of a four-terminal network can 
be expressed in terms of its scattering coefficients 
Sy, Sw, and S», and the voltage reflection coefficients 
lr, and Tr, of the generator and the load. 

Referring to figures 1 and 2, the insertion loss is 


2 
1 


P, I 
L=10 logio P,P, (10) 


:" 1 
=10 logy ;,° 
- P, ” 
where 7 is the efficiency of the network. If A and 
B denote the incident and reflected voltage waves 
shown in figure 2. The efficiency is ° 





iin P, ae E,(— I.) a” b, 2_ A, . B, "eo Py - (11) 
” 7 | ae E,- 1, A, 2_ B, a A, l — Tr, . 
The scattering equations of the network are 
B, Sy Ay “T Sie) 
; (12) 
B, SrA; T SxA,) 
Z6 - 
ry 
e 7 rs fi z re fe 2¢* Zo P "i z+ Zo 
” zy eS (i-GI+h) 
re *%5+2, "2 U-Gh) 


E e (t-GMI-G) 
k* 30" Wp eRe 











2 I-G a 2 
e 6 
R* Eri, = 425 I-@R { 1 In } 


Fieure 1. Load connected directly to generator. 


' Note that expressions of the form E+] represent the dot, or scalar, product 
of two vector quantities 
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2 2 2 
Ba} -|A2| _— [8 | 2 
P, = E2*(-I2) = ™ = _@ {11m } 


Fiat RE 2 Network inserted between gene rator and load, 


) where TP, ~p,/a, and T, ~ ayn But | Since : 
2/By S30, 
| Si. _. a (] 
B, Sho 13 ad l Sool", ; 
——9 (13) 
A 1 — Sol ; ; 
therefore. the insertion loss may also be written 
) S12 a Ty i4 R S R 
v (S— Sela] 1-187 L=20 log, 1 — Sto CU — Sal) — Siecle! (1g) 


Si.(1—P Py) 


WasHINGTON, August 28, 1953. 
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Heat Capacity, Heats of Transitions, Fusion, and 
Vaporization, and Vapor Pressure of 


Octafluorocyclobutane' 


George T. Furukawa, Robert E. McCoskey, and Martin L. Reilly 


The heat capacity of octafluorocyclobutane of 99.98,-mole-percent purity was deter- 


mined in an adiabatic calorimeter from 17 
abs j mole! 


value of 23721 
capacity peaks at 


20 abs j mole~! 

141.5 0.2°, 174.6 
slope abruptly, is attributed 
from 177° to 274° K 
to 274° K can be represented Within 


1315.906 


logyp(mm Hg) = 6.70267 r 


The entropy of octafluoroeyclobutane in the ideal gas state at 261.25 
puted from the experimental data to be 380.52 


leg~' mole~! (1 cal= 4.1840 abs j 


0.2°, 


to 270° K 
at the triple-point temperature of 232.96 
zation experiments at the vapor pressure of 590.9 mm Hg and 261.25 


+8. 778482 


The heat of fusion is 2768.2 ~ 2.0 
0.05° K. The three heat-of-vapori 
K gave an average 


Four solid-solid transitions were observed with the heat- 
214.84 
anomalous behavior in the heat capacity at about 97 
to glass transformation 
The results of this measurement from the triple-point temperature 
0.5 mm Hg in most cases by the equation 


0.05°, and 216.99 +0.05° K. The 
K, where the heat capacity changes its 
The vapor pressure was measured 


10-°T — 1.739691 * 10-57 


K and | atm was com- 
0.71 abs j deg”! mole”! or 90.95 O.17 cal 


The entropy computed from the spectroscopic assign- 


ments by Claassen and the molecular constants based on the Dy, model amounted to 87.83 


cal deg! mole 


1. Introduction 


Octafluorocyclobutane can be prepared by pyroly- 
sis of polytetrafluoroethylene * and by cyclic dimeri- 
zation of tetrafluoroethylene.’ The latter material 
has come into importance in the recent years as the 
monomer in the production of polytetrafluoroethy- 
lene plastic (Teflon), which has a remarkably high 
chemical and thermal stability. In order to obtain 
a wider thermodynamic knowledge regarding the 
fluorocarbons, fluorocarbon polymers, and monomers 
prototype molecules, calorimetric and vapor- 
pressure measurements have been made with octa- 
fluorocyelobutane. This material is interesting also 
from the viewpoint of its structure; some *° of the 
investigators have interpreted its vibrational spectra 
in terms of a planar ring model (Dy). More recent 
spectral © and electron diffraction ’ studies seem to 
favor the nonplanar ring structure (V4). There is at 
present an investigation being conducted at the 
Bureau to make spectral assignments in light of the 
evidence for a Vg structure. 

Calorimetric studies with the tetrafluoroethylene * 
and the polytetrafluoroethylene * have been re- 
ported recently by the National Bureau of Standards. 
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2. Apparatus and Method 


The heat-capacity, the heat-of-fusion, and the 
vapor-pressure measurements with the octafluorocy- 
clobutane were carried out in the adiabatic calorim- 
eter described previously by R. B. Scott, et al."° 
For the vapor-pressure measurements the calorim- 
eter valve was opened to a mercury manometer, 
which was read by means of a mirror-backed glass 


scale. The calorimeter served as a thermostated 
container. No provision was made to stir the 
sample. The pressure readings were converted to 


standard mm Hg (g=980.665 em tempera- 
ture=0° C.) on the basis that the local gravity is 
980.076 cm sec 

The heat-of-vaporization measurements were made 
in another adiabatic calorimeter similar in design to 
those described by Osborne and Ginnings " and by 
Aston, et al.” This calorimeter contained a throttle 
valve within the space enclosed by the adiabatic 
shield. During the vaporization experiments, the 
vapor was removed isothermally by controlling the 
valve and introducing a known electric energy con- 
tinuously. A nickel resistance thermometer wound 
on the tube between the valve and the calorimeter 
container served to sense the temperature of the 
vapor removed. When a desired quantity of the 
material was collected, the electric power was cut 
off and the valve closed. The heat of vaporization 
per mole, Z, was computed from the total energy 


sec 
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nput, Q, and the mass, of material collected by 


the relation 


m, 


u 
(1-3) 
u 


in which w and are the specific volumes of the 
liquid and vapor, respectively, at the vaporization 
temperature and MV is the molecular weight. 

All temperature measurements are in accordance 
with the International Temperature Seale and 
below 90° K in accordance with a provisional scale 
based on a set of platinum resistance thermometers 
calibrated against a helium-gas thermometer. 


MQ 


m 


L (1) 


3. Material 


The octafluorocyclobutane used in this investiga- 
tion was kindly furnished by the E. I. du Pont de 
Nemours & Co. Although the sample was con- 
sidered quite pure, the possible impurities were tetra- 
fluoroethylene (C,F,) and _ chlorodifluoromethane 
(CHCIF,). The material was purified further by 
fractional crystallization. 
to that described by Schwab and Wichers." The 
purity of the material so treated was determined 
prior to the heat-capacity experiments from the 
equilibrium melting temperatures and was found 
to be 99.98, mole percent. The equilibrium tempera- 
tures and the corresponding reciprocal of the fraction 
melted, F, are given in table 1. The triple-point 
temperature (obtained by extrapolating the tem- 
perature versus 1//' plot to 1//'=0) was found to be 
232.96+0.05° K. (In this paper, unless defined 
otherwise, the figure following the plus or minus is 
the authors’ best estimate of the uncertainty.) A 
second purity measurement at the end of the heat- 
capacity studies gave essentially the same result. 


Taste 1 Equilibrium melting temperatures of octafluoro- 
cyclobutane 
Mole fraction impurity =0.00613 4 T, 


K = 273.164 °C 


Reciprocal of 
fraction 
melted, 1/F 


| K 
* 232. 7587 
232. 8012 

232. 8215 

232. SA7H 
232.9114 

232. 9278 

232. 420 

» 232. On80 


Triple-point temperature, 
232.9520.05° K. Purity, 
98.9% mole percent 


* The temperatures given are accurate to +0.01° K. Wherever temperatures 
are given to the fourth decimal, the last two figures are significant only insofar 
as small temperature differences are concerned 

» Fxtrapolated 
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The procedure was similar | 








yurity determinations the system was as- 
have ideally and to form no solid solution 


In these 
sumed to 


4. Heat Capacity 


The heat-capacity measurements were carried out 
from 17° to 270° K. The observed values are given 
in table 2 and are also plotted in figure 1 to show the 
four solid-solid transitions observed. The tempera- 
tures given in the table are the mean temperatures of 
heating intervals. The temperature rise associated 
with each heating interval is also given. The cor- 
responding heat-capacity values, which are under 
conditions of saturation have been corrected for 
curvature (see footnote 16) and vapor saturation 
wherever applicable. 

The heat capacity in the temperature range of each 
transition was carefully investigated. Of the four 
solid-solid transitions, two were found very close 
together, similar to those in hydrogen bromide." 
Equilibrium was attained within these two transi- 
tions (214.84° and 216.99° K) in the usual time of 6 
to 10 min. On the other hand, the relaxation time 
for the transitions at the two lower temperatures 
(141.3° and 174.6° K) was found to be extremely 
long; small but definite temperature drifts were 
observed even 2 hr after the end of heating. On 
account of this long relaxation time, the heat capacity 
in these transition ranges was dependent somewhat 
upon the thermal history or the rate of cooling prior 
to the experiments. Thus, the results in the transi- 
tion region are scattered, up to as much as 0.3 
percent. The peaks of the heat capacities within 
the four transitions were found to be 141.3+0.2°, 
174.6+0.2°, 214.84+0.05°, and 216.99+0.05° K. 
The heat of transition was not determined for any 
given temperature; instead, the enthalpy change was 
determined for an interval of convenient tempera- 
tures, which included the transition. These results 
are given in section 5. 

Figure 1 shows an anomalous behavior in the 
region of 97° K, in which the heat capacity under 
goes a rather abrupt change in slope. Also, a longer 
time (about 30 min) was required for temperature 
equilibrium. A magnified plot of this temperature 
range is shown in figure 2. The cooling rate prior 
to the heat-capacity measurements had a decided 
effect upon the results in this temperature range. 
Slow cooling yielded significantly higher values. As 
pointed out in section 5, the increase in the enthalpy 
change corresponding to the higher heat-capacity 
values in this temperature range seems to be offset 
by an almost psc enthalpy change decrease in the 
interval 130° to 150° K 

Although not as pronounced, the heat-capacity 
behavior of the substance near 97° K is similar to the 
glass transformation found in organic glasses," high 
polymers,” and other nonerystalline solids. Similar 


" H. J. Hoge, J. Research N BS 36, 111 (1946) RP 1693. 
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TABLE 2 Observed molal heat capacities of o« tafluorocyclobutane 


Molecular weight 200.4; °K —273.16 +°% 























‘ \7 7 Cc. 47 ] ‘ \7 
RUN 1 RUN ¢ KUN ts 
h he j deg mole A kK als j deg ~' mole K the j de nole A 
7 ISA 15 4, 0158 181, SS 159. 24 7. 3718 70H 2. iat 
27 140. 01 3. GS45 189. 74 166. 45 8. O46 ss » OO 
23 191.4% 9511 198. 52 175. 39 &. 6249 “) % 2 4518 
25. 40 184.12 409 2 301 
Fusion 208. 87 189. 66 1. 774 sO. IS 2. 5287 
210. 70 193.13 1, SANS) 
4.7 199, 00 1S 211. 43 193. 57 1. 6146 
O37 8 199. 87 3. 1335 212. 83 190. 60 1. 1887 RUN 13 
241.00 200, 70 5. 1104 213. 82 202. 47 0. 7842 
246. 04 Ar. Js fi, OASe 
2h3. 35 24. 48 7. (44 lransition 22%. 41 O14 1. W004 
a) OS 207. 05 7. 5822 229. 92 92.72 14a 
Oi. 52 a. 77 34 216. 32 215.00 Pra) 
Lusion 
rransition 
RUN 2 S37. 23 ww). 62 4. 7629 
“41.97 WN. 47 4.7308 
RUN 7 
20.2 191.77 1. 6640 
211.9 195, 4 1. 7359 RON 14 
=i) 252. 55 1. a8 55. 86 SUSY 
61.45 2s4y 
Pransitior tH. 52 4.8475 SH. MH) 82.37 2. 5249 
71.20 4. 5216 NS. 62 M4. 4 2. 6105 
i 7s 273. 48 1. SSN 75. 59 4.2 91.13 sth, OF 2.4014 
; : sO. 13 4 93. 89 8s. 08 $. 1275 
Transition S483 4.5 06. H1 &U. 48 2. 2082 
sv. 20 ‘ 99. 25 92. 32 2. O47 
19.71 187. SS 2 os. 03. 57 4 101.83 On. ™) 2.1792 
7) in 18S. 5S 2 3493 WA 17 5.0130 104. 44 WS 17 2 8385 
2), til 190. 70 2. 3288 103. 07 4. 7848 107.14 100, 97 2.7717 
2&2 192. 53 4. 1002 
RUN 8&8 RUN 1 
RUN 
18. 13 17.18 0.1117 87.34 83.47 2. 1582 
Mi. 7 125. 18 9. 1535 1s. 56 18. 14 7424 SY. 65 85. 24 2. 4628 
3 14.47 143. 70 fi 19. 49 19. 1. 1328 92.09 87.05 2. 4078 
x 174.84 2 21.69 2.5 +. 2563 04.47 SA. Nt) 2. 3574 
24.04 2.17 2. (484 ¥6. SO) WO. 86 2. 3008 
rransition 27.11 11.0 2. 2921 8. 09 43.14 2. 2602 
20. 92 $5. 06 }, 3206 101. 69 8.06 2. 78 
44.24 157. 16 4. WS 2. 80 tw’. 40 2. 6231 104. 65 us OS 2. 890 
is. 2 143.83 4. 0824 6. 49 41.97 4. 5760 
4. 145. 00 5. 0648 40. 89 46. 00 4. 2323 
46.20 i). 27 6. 3786 RUN 16 
42. 43 5h. 34 6. 0745 
RUN 4 58. 16 50. 72 5. 3968 
63. 31 63. 86 4. 8057 S2. 82 79. 59 2.7401 
85. 5 82. 23 3. 5220 
111.38 105. 6. 6270 8Y. 41 S08 t. S008 
118. 42 113. 37 7. 4492 RUN 9 92. 76 87. 58 3. 2920 
126. 77 124.99 0. 2458 06. 00 w.14 4. 1907 
1 a 47.21 7.273 Wy. 16 43. 1s 4. 1057 
139. a0 210. 2¢ 2. 485 17. 45 16. 21 1. 4819 102, 22 96. 70 $. O126 
18. 83 18. 69 1. 2747 106, 89 101. O68 2. 8044 
rransition 20. 03 20. 84 1. 1343 110, 69 14. O4 4. 7038 
21.24 22. 87 1. 2789 115, 31 100. 6s 4. 5350 
142. 02 24. 1. O815 22. 57 24.4% 1. 3803 120. 30 115.56 5. 4474 
143, 2h 169. 51 1. 4173 126. 13 123. 93 6. 2149 
144.0 151.79 2. 0136 
147. 02 145. 33 2. O58S RUN 10 rransition 
149. 0% 143 2. 0686 
5 5 142.75 2. 0605 153. 42 144.17 5. G25 
146. 69 10. 1044 82. 20 79. 04 5. 7036 
159. 0% 0. @510 wH. 41 a2. 6 2.7175 
1). 67 1. 1034 ay. OY 4.79 2.445 
TSO. OS 0. SA80 01.71 aH OY 2. S828 
RUN 5 RUN Il 
151.63 @. R138 “29 87. 49 2. 9379 
157. 59 2. 81 U5. Ww 89. 05 2. 4637 
r 159. 45 2. 3639 Ww. 45 W). 72 2. 4207 
73.48 210. 46 0.9319 100, 82 4.44 2. 3531 
103. 53 U7. 38 1 OM 


Transition 


7T5.M4 } 1M. 63 
176. 71 155. 2 
178. 06 1M. 72 
176. & 158. 21 
185. 59 162. 04 
104.74 171. 42 
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behavior was reported by Huffman, et al.*' for 
cyclopentene and other compounds. 

The heat-capacity values given in table 2 are 
probably accurate to 0.2 percent, except in the 
transition regions of 141.3° and 174.6° K, the anomal- 
ous region of 97° K, and the temperature range 
below 50° K, where larger tolerances must be allowed. 


5. Heats of Fusion and Transitions 


The heat of fusion was determined in the usual 
manner by introducing a known electric energy from 
a temperature just below the triple point to just 
above it and correcting for the premelting and heat 
capacity. The results of three measurements are 
summarized in table 3. The values given in the 
fourth column are the total heats required to melt 
141.9302 g of the sample in the calorimeter. Con- 
sidering the precision of the measurements and 
various sources of error, the uncertainty of the heat 
of fusion is believed to be +2.0 abs j mole™'. 


Tasie 3. Molal heat of fusion of octafluorocyclobutane 
Molecular weight 200.04; mass of sample— 141.9302 g; triple-point temperature = 
232.06° K; °K =273.16+°C 


Heat capac 

ity and 
premelting 4H AH; 
correction 


Temperature interval Heat input 


A ahs ahs j ahs ahs } mole 

227.6435 to 237.3368 4190.9 2236.9 1. 0 274.7 
227.7655 to T34.7512 37.3 1602. 2 1965. 1 2768.8 
227.4253 to 234.8100 wi57. 7 1683. 9 1063.8 2768.0 
Mean heat of fusion 2768. 2 
Standard deviation +0.7 


H. M. Huffman, M. Eaton, and G. D. Oliver, J. Am. Chem, Sec, 78, 2911 
(1068 
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TEMPERATURE , °K 
Fiaure 2. Observed heat capacities of octafluorocyclobutane 
from 85° to 113° kh. 


TaBLe 4. Molal enthalpy changes over even temperature inter- 
vals, including the solid-solid transitions and the anomalous 
region 

Molecular weight = 200.04; °K =273.16+°C 


. ; Totalenergy Heat capacity 
remperature interval input pct ro AH 
His «— Hise « 
K abs j mole abs j mole abs j mole 
122.1832 to 150.2725 * 5095. 7 — 1010.0 4085. 7 
133.2603 to 150.6635 * 3739.7 +346.0 4085. 7 
Mean 4085. 7 
129.2400 to 151.6058 © « 4405.2 328.8 4076.4 
Fy x — Hire « 
170.9387 to 180.8963 * Sor. 3 +H.0 SOR. 3 
168.7525 to 177.2341 * 4768. 1 +237.8 5005. 9 
Mean 5007.1 
Hone «— Hoe « 
209.4137 to 220.8877 * 4878.4 -278.9 4500.5 
212.4227 to 222.1278 * 4529.9 +70. 0 4599.9 
Mean 4599. 7 


Hye: « — Hee « 


87.1084 to 105.4621 * 1659. 0 — 289.9 1368. 1 
89.9824 to 105.7577 * 1450.4 —81.2 1369.2 
Mean 1360. 2 
88.4184 to 106.1372 * « 1626. 9 —247.7 1379. 2 
87.7140 to 103.7215 « 1446.5 —67.3 1379. 2 
Mean 1379.2 


* Heated in increments over this temperature interval under adiabatic con- 
ditions 

>» Heated continuously over this temperature interval under adiabatic cor 
ditions. 

¢ The sample was cooled slowly over 1 week from 230° to 8° K prior to the 
heat-capacity Measurements 
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\s mentioned in section 4, in the temperature 
nee of the four solid-solid transitions, small incre- 
nts of heat were added to obtain the shape of the 
at-capacity curve. Instead of attempting to 
‘ermine the heat of transition at a given tempera- 
ire. the enthalpy change was determined over an 
nterval of convenient temperatures. This was done 
both by summing the incremeats of heats added in 
the determination of the shape of the heat-capacity 
and by heating continuously through this 
The results of both experimental methods 
tre summarized in table 4. The values given in the 
hird column are corrections used to reduce the 
nthalpy changes to the even temperature intervals 
idieated. The anomalous region was treated in a 
similar manner by summing the heats added in the 
nterval from 90° to 105° K. In one series of meas- 
irements (see table 2, runs 15 and 16), the material 
vas cooled slowly over a period of 1 week from 230° 
to 80° K. The enthalpy change, as well as the heat 
capacity in the interval 90° to 105° K, was higher. 
The enthalpy change from 90° to 105° K was 10.0 
abs j mole~' higher. Also, there was found a decrease 
in the enthalpy change in the interval 130° to 150° K 
of 9.3 abs } mole! 


irve 
! terval 


6. Heat of Vaporization 


e 

Three heat of vaporization experiments were 
made at 261.25° K (p=590.9 mm Hg). The ex- 
perimentally observed quantity, \/Q/m, the con- 
version factor, 1—wu/u’, and the heat of vaporization, 
L, are given in table 5. The specific volumes, u 
and u’, for this material were kindly furnished by 
the E. 1. du Pont de Nemours & Co. Considering 
the precision obtained and various sources of error, 
the uncertainty of the heat-of-vaporization value 
s believed to be +20 abs j mole 


laste 5. Molal heat of vaporization of octafluorocyclobutane 
at 261.25° k 
I sure = 590.9 mm Hg; molecular weight = 200.04; °K =273.1¢ ( 

WQim 1-1 I 

aha j mole ahs mole 
23, 842 0. 98550 23, 735 
23, 823 ws) 23, 716 
23, 818 W550 23.711 
Mean 23, 721 
Standard deviation +6 


7. Vapor Pressure 


\s there was no provision made to stir the sample, 
several series of vapor-pressure Measurements were 
made under different conditions io make certain 
that temperature equilibrium was attained. The 
neasurements were made going down the tempera- 
ture seale as well as up. The close agreement in 


the results given in table 6 indicates that good 
emperature equilibrium existed. 
he 


The results above 


triple-point temperature (232.96° K) were 


| 
| 
] 
| 


15 


fitted to an equation given by 


oe ~ 1315.906 
login p(mm Hg)= 6.70267 — r - 


8.778482 * 10° 7 — 1.739691 & 10°T" 
The observations are compared with this equation 
in table 6. The equation gives the normal boiling 
temperature as 267.17+0.01° K and the tripie-point 
pressure as 142.8 +0.3 mm Hg. 





TABLE 6 Vapor pressure of o¢ tafluorot yclobutane 
° K =273.1¢ ( 
7 p p p p 1 p p p I 
Series I K mm mm mm He 
He He 
241.32 226.3 226.5 2 
K mm mm mm He 246.90 303.3 13, 3 “0 
b _He He 253.01 | 406.5 406.5 i) 
267.89 782.6 782.9 0.3 250. 51 47.3 | 547.6 
270. 95 SSH. 2 BADD 23.88 661.9 H2Z2 
273.63 Q84.8 G84.2 t 206.65 | 743.7 743.9 4 
Series IT Series III 
179. i¢ 177. 36 1¥ 
191. 04 191.82 8.2 
Di. 23 1 106. 86 wY¥ 
215. 30 4: 197.17 11.0 
223. 0 77 26. 70 21.7 
230. 81 125 213.09 5.9 
232. 91 142.4 142.5 0.1 221. 67.7 
233.16 144.6 144.5 ! 271.94 921.4 921.3 ot 
8. Experimental Entropy 
The experimental d : 
1¢ experimental data were used to calculate the 


entropy of the octafluorocyclobutane in the ideal gas 
state at 261.25° K and 1 atm. The results of the 
various calculations are summarized in table 7. In 
the temperature intervals in which the transitions 
occurred, the entropy change was evaluated by sum- 
ing the various experimentally observed CAT/7’s. 
where C is the observed heat capacity, AT the temp- 
erature rise, and 7’ the midtemperature of the heat- 
ing interval. The entropy change so obtained was 
corrected to the even temperatures indicated in 
table 7. Below 180° K, the results from the slow- 
cooling experiments were used to evaluate the en- 


TABLE 7 the 


Summary of experimental molal ¢ ntropy af 
or tafluoroc yclobutane 
Molecular weight = 200.04; °K = 273.164 °C 
aba j) deg male 

Siso, Debye 1 3840.04 
AS °, numerical 7S. 2440. 16 
5 So transition 14. 15). 10 
AS numerical 24. 03-40. 05 
AS , transition 2. 1440.06 
AS °, numerical 15. 67 20.04 
4S eo®, transition 75. 65440). 08 
ASiwe ©, numerical 2h. 4-4-4). 05 
AS; ~»°, transition 21.3240. 02 
ASom «°, numerical 0. #240 02 
ASng.eee, fusion, 2768.2/242.06 11. S40 O1 
ANSing.® , numerical 23.2240 4 
AS ©, vaporization, 23721/261.25 0) 80-40 O48 
AS , gas imperfection 0. 67 

AS ©, compression 20 

Entropy, ideal gas at 261.25° K and 1 atm oe0) 0.71 











tropy. The gas-imperfection correction was com- 
puted on the assumption that the octafluorocyclo- 
butane vapor can be represented by the Berthelot 
equation of state. The critical constants used were 
those reported by Renfrew and Lewis.” 

The uncertainties given in table 7 were based on 
the precision, the probable nonreproducibility of the 
physical state below 180° K, and various other 
sources of error. The absolute entropy obtained is 
probably in error by more than is given in the table 
because the results of the measurements indicate a 
certain degree of glassy state or randomness in the 
solid at the lower extreme of our measurements. 


9. Entropy From Molecular and Spectro- 
scopic Data 


The entropy of the gaseous octafluorocyclobutane 
was computed from the molecular and spectroscopic 
data based on the Dg, structure. The results are 
summarized and compared with the calorimetric 
entropy in table 8. The rotational entropy was 
computed directly from the principal moments of 
inertia given by Smith, et al.” in which the C—C 
and C—F bond distances and C—C—C and F—C— F 
angles were 1.53 A and 1.36 A and 90 deg. and 112 

2M. M. Renfrew and E. E. Lewis, Ind. Eng. Chem. 38, 870 (1946 


» DD. C. Smith, J. R. Nielsen, L. H. Berryman, H. H. Claassen, and R. L 
Hudson, NRL Report 3567 (September 15, 1949 





| 
| 
| 
| 
} 





The 


wer® 
[,=I,=796 X10-" g-em? and /,=95110-" g-em?. 
The vibrational entropy was calculated on the basis 
of the assignments by Classen (see footnote 5). 

On the basis of the existence of residual entropy 
in the calorimetric work, the discrepancy in the two 


deg., respectively. moments of inertia 


results given ia table 8 would be increased. Work 
is now in progress at the Bureau for the reassignment 
of the frequencies on the basis of the recent evidence 
for a V, structure. 


Molal entropy of octafluorocyclobutane at 261.25° K 
from molecular and spectroscopic data 


TABLE 8, 


Molecular weight = 200.04; °K =273.164+°C; 1 cal =4.1840 abs 


cal deq™' mole 


S, translation 41.13 
S, rotation 25. 65 
S, vibration 21. 05 
SN, total, ideal gas at 261.25° K and | atm 87.83 
S, calorimetric 0. 95 


The authors are indebted to E. I. Du Pont de 
Nemours & Co. for the sample of purified octafluoro- 
cyclobutane. 


WasHINGTON, July 24, 1953. 
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A Characterization of Normal Matrices’ 


Alan J. Hoffman and Olga Taussky 


A matrix A is called normal if AA* 
matrix of A 


A*A, where A* is the transposed and conjugate 
It is known that for a pair of commuting matrices, A, 
ordering of the characteristic roots ay, » Ge 
polynomial p(A,B) has as characteristic roots the numbers p(a;.8 

general, weaker than commutativity, but does imply it if B=A*. It 


B, there exists an 
, 8, of B, such that every 
This property is, in 
is shown that this 


of A and §;, 


property already implies commutativity of A and A* if it is assumed to hold for only one 


polynomial, provided the latter is suitably chosen 


are examined for their suitability 


Normal matrices can be characterized in many 
avs. In the original definition [1]* a matrix A 
called normal if AA*=A*A, where A* is the 
ransposed and conjugate matrix of A. It is known 
that for a pair of commuting matrices A, B, there 
cists an ordering of the characteristic roots a, . . ., 


of A and 86, 


Ay 
. ., B, of B, such that every poly- 
nomial p(A,B,) has as characteristic roots the 
numbers p(a;8,), see [2]. This property in 
veneral, weaker than commutativity, but does imply 
it if B=A*. This follows from the fact that two 
matrices with this property can be transformed to 
upper triangular form simultaneously by a unitary 
similarity transformation [3]. Even if the property 
s assumed only for linear polynomials p(A,A*), 
commutativity follows. This is a consequence of 
the fact that A is the sum of two normal matrices, 


Is, 


' A+A* , A—A* 

; a 
For the matrices (A+ A*)/2, (A—A*)/2 will again 
have the above property with respect to linear 


polynomials; this, however, implies that they com- 
mute [4, 5]. Henee A and A* commute. 

It will now be shown that the above property 
already implies the commutativity of A and A* if 
it is assumed to hold for only one polynomial, 
provided the latter is suitably chosen. Two special 
results in this direction were obtained earlier [6]; 
the polynomials considered then were A+ <A* and 
AA*, and it was assumed that the special ordering 
of the characteristic roots implied that 8,—a,. 

Certain polynomials of first and second degree are 
examined for their suitability. The result is that 
any polynomial of first degree is suitable, while not 
every polynomial of second degree is, although a 
large class among them is suitable. 

In what follows we assume that A=(a,) is an 
nn matrix with complex numbers as elements. 

THEOREM 1. Tf a, _ a, are the eigenvalues of 

nn matrix A, a and B are two complex numbers 
lifferent from 0, and the eigenvalues of aA+BA* are 


} 


| 
¢ preparation of this paper was sponsored (in part) by the Office of Scientific | 
earch, United States Air Force | 


| 
Figures in brackets indicate the literature references at the end 


paper 
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Polynomials of first and second degree 


+-Bap;, where P is a permutation of the integers 
1,. . ., n, then A is normal 

Proor. It is clearly no restriction to assume 
8=1. Further, since the hypothesis and conclusion 
of the theorem are concerned with eigenvalues of an 
operator A in unitary n-space, we may assume that 
an orthonormal basis has been chosen in the unitary 
n-space, so that the matrix A is (upper) triangular; 
i. @., @.=0 for i >k, and ay—a;. We shall show 
that 


aa 


a.=—0 for i+k. (1) 
which implies that A is normal. 
The argument is based on the fact that the second 


elementary symmetric function S, of the eigenvalues 


V1, Yn Of aA+A* may be computed in two 
wavs. First, since 
2 aa; + ap (gaz I, n) 
we have 
n 
P ’ ‘ Vy _- = - = ‘ 
28S.=3;5 | a acy ay, + A py py + A pA py + Apyay)|. (2) 
j 1 ki 


Second, since S, is also the sum of all principal 
22 minors J, of aA+A*, we have 


n ” 
. ’ ‘ ‘ ‘ 2 -- 
2S.=>>5 S* Ma=>°* > {ac xj ary, + Qe, 
i=l kei 1 kei 
+ af aa, a;,a,—a 1 y)]. (3) 


As i and & range over all possible distinct ordered 
pairs of elements of the set {1, . . ..»},so do P, and 


P,. Hence if we subtract (3) from (2), we obtain 
" 
> > > > (Qj py > A py ay aya, a; +A,Ae)—0. (4) 
i=l kxi 
P , 
But @y¢”>0. Therefore, if we show that 
n 
RS S* (aya@pe + appa, — aya, — a;a,) >0, (5) 
i=l ket 


then (4) will imply (1). 














Write Since a,44,,>0, the only way for (12) and (13) to 
ay Aj (gm ], . «. <8), be consistent is for 


where A, and u, are real 4,4, —0 for all i +k; 
Then (5) becomes 
i. e., A is normal. 


" -¢ se esse _ ss ’ arg . , 
25) SIA a oe en (6) a use essentially the same argument t 
THEOREM 3. Let a, .. ., @n be the eige nevalues of th 
We shall prove an nXn matrix A, P a permutation of the integer 
l, . n, and ay, ay, by, by, e), 2 compler numbers 
SS A ..)>0 (7) Let the ergenvalues of 
—_ _ Pi k _ ° ‘ 


aA T a,A* T b, A? T 6.4* T ¢,A A* T C2A* A (14) 


The same argument will apply to the u’s, and P 
henee (6) will be established. Now for fixed 7 ” 


SN Npe— Ag) = Ag DO (Ape — Aa) = AA — Api) 00; + 2 py + Dyce? + boa}, + Cyaap 
“T C2ap a, (7 aie ee Gee 
Hence (7) is equivalent to If, in addition 
¢, +e.#0, (15) 
Sy (.— dv.) >0 (8) the nA is normal. ; 
es What if condition (15) is omitted; i. e., if ¢,=e,” 


Then theorem 3 may not hold. There are two Bi 
Inequality (8) is a special case of a known in- | C88¢s to consider: 
equality ({7], theorem 368), which states that if 
a@,>.. .>@,, 6,> .>b,, and P is any permutation 
of {1,. . . ,n}, then 


(1) ¢,=0, (i) ¢, #0. an 


In case (i), if 6b; =b2=0, then the theorem reduces to 


6 
py? SSR O (9) | theorem 1, so A is normal (unless a,a.=0, where the 
— a) Pi- . . nw . 
theorem fails). The theorem fails unless 6,=)b.—0 
Set a,—6,—X,, and (8) is an immediate consequence A counterexample is furnished by 
of (9) | by 
Turorem 2 Tf a, . a, are the eigenvalues of e d Wa 
an nxn matrix A, and the eigenvalues of AA* are ne 
Oya p,, , AyGp,, where P is a permutation of the 0 0 
integers 1, . . ., n, then A is normal. : 
Proor. As in theorem 1, we may assume A is where dO is arbitrary, ¢ -a,/b, (if 6,40), or 
il lar We shall a *; c —@2 b, (if b, #0). 
upper triangular e shall compute trace AA* in F ef seueti ' ' “ee 
~ ne . mm is some —_ 
two ways. First, since trace AA* is the sum ot the oli aed @, fin “ fa S em I | ve 
. 7 . . . ’ , ey . e . Te “rea 
eigenvalues of AA*, we have by hypothesis Cele age a Se ie Ce ee ee h 
= tive, the theorem fails, as one can see from the 


counterexample 


n 
trace AA* aap, aap, > 0. (10) 0 a 
i=l | 
ii | 0 o/, 
Second, since trace AA®* is the sum of the elements | where a’¢@——a,a 
. j * "e eo . - . . ene 
on the diagonal of AA*, we have On the other hand, if (14) is formally Hermitian, 
. " the theorem does hold. Specifically, we have 
trace AA*= Slay t+ 3) Sandu. (11) THEOREM 4. Let a, . . .,a be the eigenvalues of 
i=l i=) kei an nXn matriz A; by P be a permutation of the 
integers 1, . n; a, 6 complex numbers, c, d real 


Write a,;= pe’, where p;>0. We have, from (11) 


and (10) numbers not both 0; let the euge nvalues of 


n Nn n i! a a 8 we hae . 
Dei + Dy Didudu=D pen. gy | APSE ee tee ee 
i=l t=1 kei i=l 
+dA*A (ec, d real, not both 0) (16) 
But from inequality (8), we see that be 
Dota Dover. (13) | Geet Gains + bad + bap, + cara, 
i=l 


+ dap i + + eee 
This theorem can also be proved from known inequalities, see [8], between ° . la iy ( l , 
the eigenvalues of a matrix A and the eigenvalues of AA* then A is normal. 











oor. If e#—d, the situation is covered by 
em 3. Ife d, it is clearly no loss of gen- 
y to take c=1, d l. Now order the 
ivalues of A, so that 

P (aay, 1 bai) 4 > PR iaa, + be ) (17) 
We may assume that A is upper triangular, and 
at the eigenvalues of A appear on the main diagonal 


the order a, @,. . ., @p. 
Then let B stand for the Hermitian matrix (16), 
d denote its eigenvalues by 8,> . . > B, 


It has been shown ( [9], theorem 1) that 
5,,<68 


bi, +b <8, 4-8: 


1S) 
4 h . — p T B, 
But 
Rix 2P (aa, +- ba?) Sa, Su, (19 
and 
5 max|/? (aa dap be bad }<2?P (aa, 1 ba), 
(20 
by (17). Since ay,.4,>0 (k=2 , n), the only 


wav for (19) and (20) to be consistent with the first 
inequality of (18) is if 


a 0 (-k—2, . i) (21 


Now, in view of (21), it follows that 


by, bo» 2Pla(a, tT ay) b(aé T as) | T > \42.4cx, (22 





and 
> =e = 
p Bp» max [?! a(a,+a + Aa\ap ap 
’ 
b(n 2 h( 2 at 99 
0\ a aw) O.aG y ay =>) 
Op , , 
- 22a a + ay) +-b( a? ay)! 


The only way for (22) and (23) to be consistent 
with the second inequality of (18) isifa,,—0 (k=3, 

o «6 gts 

Continuing in this way, one sees easily that (+k 
implies a,.—0; hence, A is normal 


The authors thank J. Todd and H. Wielandt for 
helpful comments 
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Prediction of the Likelihood of Interference at Frequencies 
of 30 to 42 Megacycles in Alaska 


T. N. Gautier, Jr., and C. J. Sargent 


The likelihood of interference with 


very-high-frequency 


networks in Alaska from 


stations operating on similar frequencies in other parts of the world, as a function of season, 


sunspot number, and time of day, is presented 


computations for other communication links. 


This example serves as a model for similar 
Calculations are based upon regular reflec- 


tion from the F2 layer, and data are supplied for estimating the probability of interference 


from sporadic-F reflections 


1. Introduction 


\t the request of the Civil Aeronautics Administra- 
ma study was made of the likelihood of interference 
th the operation of a proposed very-high-frequency 
ne-of-sight network in Alaska caused by ionospheric 
ropagation of signals between stations within the 
vetwork and of signals from distant stations outside 
the network. The study was made for operating 
frequencies of 30, 36, and 42 Mc. The likelihood of 
interference at these frequencies is inferred from the 
likelihood of occurrence of maximum usable fre- 
quencies (MUF) equal to or greater than these 
frequencies calculated for propagation paths between 
stations within the network and between distant 
points and two representative points in the network. 

This specific study is an example of a communica- 
tion problem of general interest, and is presented in 
full to serve as a guide in solving similar problems. 

Nome, Alaska, and Annette Island, Alaska, near 
the extremes of the network, were chosen as the 
points within the network for the study of interfer- 
ence from distant stations. 

The likelihood of the MUF exceeding 30, 36, and 
12 Me was calculated for three classes of propagation 
paths: (1) F2-layer paths from stations at distances 
greater than 4,600 km from Nome and Annette 
Island, respectively. (2) F2-layer paths less than 
+,000 km in length, including paths between network 
stations. (3) Sporadic-F paths between stations of | 
the network. 


2. Basic F2-4000 MUF Prediction Charts 


Using standard procedures, F2-4000 MUF contour 
charts similar to those appearing in the CRPL-D 
series [1] ' were drawn for latitudes 30° to 70° north 
in the I-zone, for sunspot numbers 60, 80, 100, 120, 
and 140, and the months June and December. It | 
was not necessary to prepare charts for sunspot 
numbers less than 60 for the problem under con- 
sideration. 

Sunspot numbers referred to are 12-month aver- 
of relative sunspot numbers. The average 
maximum 12-month average sunspot number for the 
past 10 cyeles is 105, with a standard deviation of 24. 
‘he maximum for the present cycle, which occurred | 

1947, was 152. 


aves 


ures in brackets indicate the literature references at the end of this paper. 
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3. Interference from Distant Stations 


Using the charts referred to in section 2, figures 1 
to 34 were prepared. Figures | to 16 are for recep- 
tion at Nome, whereas figures 17 to 34 are for recep- 
tion at Annette Island 

Curves on these charts show the predicted varia- 
tion with sunspot number and time of day of the 
percentage of occurrence of conditions for reception 
at 30, 36, and 42 Me by F2-layer propagation from 
transmitting stations located at distances of 4,000 
km and over, in selected azimuth ranges measured 
clockwise from north at the receiving station. The 
azimuth ranges for paths terminating at Nome are 
53° to 117°, 172° to 254°, and 254° to 299°. The 
azimuth ranges for paths terminating at Annette 
Island are 67° to 157°, 218° to 277°, and 277 
322°. The world areas included in these azimuth 
ranges are shown in figures 35 and 36. Charts for 
summer and winter conditions are presented. 
Charts on which the predicted occurrence is less than 
10 percent for sunspot number 140 are omitted 
Equinox conditions are approximately intermediate 
between those for summer and winter. 

As an example, in figure 1 the chart for Nome, 30 
Me, June, and azimuth range 53° to 117°, conditions 
for reception of 30 Me at 0500 GCT on 90 percent 
of the days of the month on the average should occur 
for a sunspot number of 135. For sunspot number 
70, reception should occur on only 10 percent of the 
days on the average at 0500 GCT, whereas for sun- 
spot number 120, reception on 10 percent or more of 
the days should occur between about 1830 and 0730, 
with about 75-percent occurrence at 0500. 

The charts were constructed as follows: 

1. The 2,000-km circle around each 


to 


receiving 


| station shown in figures 35 and 36, and the selected 
| azimuth intervals, were traced on transparent paper. 


2. Daily intervals for reception on more than a 
specified percentage of days were found by placing 
these sheets on the predicted 4,000-MUF charts, 
and noting the earliest and latest times at which the 
appropriate MUF contour intersected the segment 
of the 2,000-km circle included in the given azimuth 
range as the transparency was moved horizontally 
from left to right across the MUF chart. 

For reception on 50 percent of the days, the 
appropriate contour is that for MUF equal to the 
given frequency. The dispersion of daily values of 














the MUF is such that contours for MUF equal to 


1/(1+0.15) and 1/(1—0.15) times the given 
queney correspond approximately to 10- and 90- 
percent reception, respectively, [2, see. 6.6, i). 
Assuming that the distribution of percentage of 
MUF is Gaussian, contours for MUF equal to 
1/(1+0.08) and 1/(1—0.08) times the given fre- 
queney would then correspond to 25- and 75-percent 
reception, respectively. 

3. Contours in figures 1 to 34 were then plotted 
by interpolation between the end points of the inter- 
vals determined in step 2. Figure 17(b) is discussed 
below 

Although the predictions are based upon an 
assumed ring of transmitting stations 4,000 km 
distant from the receiving station extending over 
the indicate] azimuth range, stations at greater 
distances in the same azimuth range are also in- 
cluded, provided the MUF at the control point 
nearest the transmitting station is greater than the 
given frequency. This will be true, usually, for 
transmitting stations in the land areas included in 
the azimuth ranges, except those west of the receiving 
station during the early hours of the predicted daily 
reception period, and those east of the receiving 
station during late hours of the predicted daily 
reception period, at distances well over 4,000 km 

Actually, portions of the 2,000-km circles around 
the receiving stations lie outside the I-zone. How- 
ever, errors introduced by using I-zone predictions 
for these segments are believed to be less than other 
uncertainties in the predictions. The error in the 
position of portions of the curves on these charts may 
be of the order of 10 or 15 sunspot numbers. 

It should be emphasized that these charts show the 
occurrence of conditions for interference from one or 
more of a number of distant stations operating on the 
same frequency and distributed over the entire 
azimuth range. The interference from a single sta- 
tion would be somewhat less than this. For example, 
the predicted interference at Annette Island from a 
station in Washington, D. C., operating at 30 Me in 
June is shown in figure 17(b). This may be com- 
pared with figure 17(a), which shows the predicted 
interference from stations throughout the azimuth 
range 67° to 157° 


4. Interference from Stations Within the Net- 
work and Other Stations at Distances 


Less Than 4,000 Km 


Figures 37 to 49 show the predicted occurrence of 
conditions for interference at 30, 36, and 42 Me, via 
F2-layer reflections from stations at distances less 
than 4,000 km as a function of distance and time of 
day for paths centered on 60° N. latitude, the approx- 
imate latitude of the middle of the network. Charts 
are given for June and December and for sunspot 
numbers from 60 to 120, in steps of 20. They were 
prepared by using the MUF values at 60° N. on the 
F2-4000 MUF charts and assuming a ratio of 4000 
MUF to zero MUF equal to 2.9, and interpolating 
by means of the nomogram on page 85 of reference 


fre- | 


TEs above 7 


[2]. Conditions for equinoctial months are inter 
mediate between those for summer and winter 
When the predicted occurrence was less than 10 per 
cent at all times of the day for distances less thar 
4,000 km, the chart was omitted. 

It may be noted that the predicted occurrence fo 
distances less than 2,000 km, the maximum distance: 
within the network, exceeds 10 percent only in winte: 
near sunspot Maximum. 

Figure 50 shows the percentage of occurrence by 
seasons of vertical incidence sporadic-F reflections 
(fEs) at frequencies above 7 Me, observed at Anchor 
age, Alaska, averaged over 4 vears, May 1949 through 
April 1953. Occurrence of conditions for propaga- 
tion of 30, 35 and 40 Me, based on the occurrence of 

7 Me may be estimated as follows: 

The logarithm of the probability 7? of occurrence of 
vertical incidence sporadic-- reflections at fre- 
quencies higher than f is roughly a linear function of 
the frequency [3]. Thus 


logiy P=a+ bf, 


or 


logio ( P, P,) 


PP: 


b(f:— fa) 
10° 


This does not apply if P is close to 1.0. 

Analysis of Anchorage {Es measurements indicates 
that the occurrence of fEs above 5 Me is approxi- 
mately 2's times that for fEs above 7 Mc. From the 
preceding relationship B=(0.4)', so that P, 
P(0.4)%~-"? and the probabilities of occurrence of 
vertical incidence Es reflections above 6, 8, 9, 10, and 
11 Me are 1.6, 0.6, 0.4, 0.25, and 0.16, respectively, 
times the probability of occurrence above 7 Me. 

Assuming régular E-laver MUF factors, ks MUF 
equal to 30, 35, and 40 Me, corresponds approxi- 
mately to fEs equal to 6,7, and 8 Me, respectively, 
for a distance of 2,000 km, to Es equal to 7, 8, and 9 
Me, respectively, for a distance of 1,500 km, and to 
fEs equal to 9, 10, and 11 Me, respectively, for a 
distance of 1,000 km. 

The approximate relative frequencies of occurrence 
of ks MUF greater than 30, 35, and 40 Me, for 
distances of 1,000, 1,500, and 2,000 km (i. e., ratios 
to the occurrence of fEs above 7 Mc), based on the 
above relations, are summarized in table 1. 

Referring again to figure 50, Anchorage fF» above 
7 Me shows little seasonal variation, which is in 
accord with observations at other stations in the 
auroral zone. Little is known about the sunspot- 
cycle variation. 


TABLE 1 


Distance 


Frequency 


1,000 km 1,500 km 2,000 km 
Me 
uw o4 LO 1.6 
a 25 0.6 1.0 
Pt) 16 4 o6 
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Fiaures 35 and 36. 


Maps showing geographical areas included in selected azimuth sectors. 
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Figures 37 to 41. Predicted skip distance for F 2-layer reflection at 60 





N. latitude, I-zone. 























IN KILOMETERS 


DISTANCE 


DISTANCE IN km 






























































































































































































































































































































































FiG. 42 FIG. 45 
30Mc SUNSPOT NO 100 JUNE 36 Mc SUNSPOT NO. 100 DECEMBER 
4000 - 4000 
i NUT ) " 
| | 
| | 10% | = Lj 
3500-——> > + + +~—+ 35 + 
TTT TT | | | 
I} | | | | | | | 
sooo aul 3000 | + t — + + 
LEE EE 
° 1 | | | 
FIG. 43 _— 
ani 50Mc ; SUNSPOT NO !20 > FIG. 46 
36 Mc SUNSPOT NO. i20 OECEMBER 
| 4000r—y T | f 
340 i + av . | | 
+ v PTTL | 
Oo / 3500 + + | | + + + 4 | 
Be }] 7] | | | | | 
3000 + J | | 
7 
Nie Y 3000 oe + 4 | 
2500 + i , + N 
| 2500 } aes + + | 
‘a rT 
i a a 
FIG. 44 FIG. 47 
30Mc SUNSPOT NO 140 JUNE 36 Mc SUNSPOT NO. i140 DECEMBER 
: \ | | \ | | ] | | l 
35 \ N 7 ] 4 3506 t \ t 4 | 
/ ' NS | ; 
3000 2 SY 3000 coe hex) + 1 
2500 + + J 500+ — ~ - 
+ TT — 
Wes / \ »/| 
— | j at, 
2000 + | + } 
is —_—— 1s00 1 — 
O 4 6 28 0 2 my * * 20 22 00 00 02 _ Os oe 10 2 ‘4 6 r 20 22 oO 
LOCAL TIME AT THE MIDPOINT OF THE PATH 
Figures 42 to 47 Predicted sh p distance for F 2-layer reflection at 60 NV. latitude, l-zone 
| 
Fig 48 Fig 49 
36 Mc SUNSPOT NO 140 JUNE 42 Mc SUNSPOT NO. i40 DECEMBER 
4000 TT TT TT TT TTT TT rT 
| | | | \| | | | | | | | | \\ 
| | | 
3500 | ai ani + T cos a 
Baan Net | | \ | | 
| | 
10 %, 
3000 + ’ 
| | | | | YS | | | 
2500 LI a8 aE FSSens | LLL 
2 O84 € 8 4 8 20 22 OO 0O 02 O04 O6 O8 10 12 i4 16 1®@ 20 22 00 
LOCAL TIME AT THE MIDPOINT OF THE PATH 
Fieurnes 48 and 49. Predicted skip distance for F2-layer reflection at 60° N. latitude, I-zone. 


30 























MAY JUNE JULY AUGUST 






































6 MARCH APRIL SEPTEMBER OCTOBER aed 
© 
Oo -—-6 ae 
- 
za 
WwW 
YU 
x 
wW 
a 
| NOVEMBER DECEMBER JANUARY FEBRUARY 
| y 
: ~~ 
] 
4 2 oom 
@) 6 12 18 24 
HOURS, LOCAL TIME 
RE OO Diurnal variation in percentage of time of occurrence of sporadic E reflections in exce 0 V lnchorag TZ 


Vay 1949-April 1953 





WaAsHIneton, July 1953 











ae 
fie 
tre 
at 
ol 
eX] 
ne 
ar 
the 
the 


a 


Tables for Use in the Interpretation of Paramagnetic 
Behavior Below 1° K; for the Chromic Alums (J=3/2) 


Ralph P. Hudson and Charles K. McLane 


of interpreting susceptibility 
experiments, | 
recounted 


this procedure 


Vv mea 


mperatures briefly 


> CarT 


chromic 
that the 
ion 


The low-temperature behavior of the 
ims 1s largely determined bv the fact 
fold degenerate crround level of the Cr 
by the crystalline electric field Stark effect 
two doublets separated by an energy The 

of Hebb and Purcell [1] ' yields expressions 
f« e susceptibility, X and the entropy, S, as 
inctions of the absolute temperature, 7, and this 
Thus one may obtain a relation between 
S and xy, for ans In practice one usually 
emplovs the “magnetic temperat ’ T* instead of 

and 7 c/x is the Curie constant per 
ubie centimeter. 

Now experiment leads to a x S 
/ versus & relation, for NS mav be calculated quite 
accurate ly simply irom the conditions of macnetic 
field and temperature at beginning of the isen- 
tropic demagnetization process, and xo is Measured 
One therefore has only to find that value 
of 6 which leads to an S versus 7* curve fitting the 
experimental points. Then over the range of agree- 
ment one may suppose the Hebb and Purcell xo T 
and S(7) formulas be individually correct, and 
the former leads at the establishment of 
the absolute temperatures 

Table 1 gives, effectively, corresponding values of 


T’é6 and 7*/6 ecaleulated from the Hebb and Purcell 
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! Figures in brackets indicate the literature references at the end of this paper 


282952—54———3 


= of the thre 


tained in adiabat 
lin order to derive 
to tables that 


tion 


tor the cas 


'/6, 
th] 1hT 
[at 8 )+(3—%, 


Tables 2 to contain the data necessary 
culating the “experimental entropy,” that 
magnetic field at the initial temperature 
field-free entropy as a function of kT 6 

Considering first the case of the entropy in zero 
field and at low temperatures, two effects contribute 
to a reduction in the entropy below the “ideal 
value of Plog, 4, namely, the Stark splitting and 
Interaction The theory treats these two 
contributions separately, and they additive 
The “Stark entropy” is calculated from the partition 
function, Z,.=—2(1+-e7**"), and given in table 2 
Table 3 function where z=kT'6), 
the quantity Q 2 defined in 
reference [1]. then 
given by 
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ye 
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cal- 
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of r 


being 


lists as a 


dida as 
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Is 


The magnetic entropy, S,, 


(2) 


da 
0.0211 deg and 0.0194 deg for the 
methylammonium alums, respec- 


S,,/R 


6 


where 7// 3c 
potassium and 
tively Hence 


S R log, 1 S R Stark 7 S R. 


log, 4 
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Table 4 lists values of the entropy, for correspond- 


ing values of H/T, calculated by means of the 


“Brillouin formula.” 


S/R),=acotha 
where a==pll 
9.27 10-** en 
1.3805 107-"° 


Bi: 


facoth4a 


For the case J=3/2, 


log, sinh 4a 


the Bohr 


one has 


log, sinha, 
3) 


magneton, 


iu, and & is the Boltzmann constant, 


erg deg 





In order to derive an accurate value of 6, however, 
one must take cognizance of the fact that 


the 


Brillouin entropy (table 4) is only appropriate to the 


ideal case of a fourfold degenerate ground level. 


applied to the entropy reduction, log,4 
relatively large. 

An entropy formula has been derived [3] that takvs 
this effect into account, but it is more unwieldy th 
the Brillouin formula (a straightforward function of 
/1/T) as it is a function of H, T, and 6. For thy 
chromic alums the value of 6// is always close to 0.25 
deg, and starting temperatures (that is, bath temper- 
atures) will, in general, be of the order of 1.1° to 
1.2° K. Beeause AS/R<S,/PR, the Brillouin an 
exact entropy curves run very close together, and 
the following approximate method of calculating 
new AS/P under slightly changed conditions of 
and/or 7 (say 6, 71), given a table of exact values 
for conditions 69, To, is found to be very useful. 

Let A=6/2, a=A/T, and B=yll/A. One may ther 
take the expressions for both the Brillouin entropy 
and the exact entropy and obtain their derivatives 
with respect to a@. Denoting derivatives by 
primed symbol, we then have 


Sp R. 


(S,/R), ~=(S,/R)q S,/ Ra. X(a—a Brillouin) (4 
and 
(S/R) = (S/R). +(S’/ Ra X (a— a (exact) 5 
Subtracting (5) from (4 for a positive quantity 
we find 

AS/R), =(AS/R),.+(AS"/R)a. X(a— am 6 
to a satisfactory degree of accuracy. It must be 


enphasized that eq (4) and (5) are only approximate 
Higher derivatives have been neglected, but thes 
effectively disappear in the subtraction, which leads 











lhe effect of the ery stalline electric field is to modify to eq (6). 
the pattern of the level splitting in a magnetic field, In table 5 the entropy values have been calculated 
and for small fields correction, AS/R, to be! for various values of AZ, and with 6/k=0.27 deg 
TABLE 4 Entropy as a function of a uH/kT) calculated from the “Brillouin fo ila”: S/R=a cotha fa coth 4a-+- log, s 
fa —log, sinh a 
lables for J=1/2, 3/2, 5/2 ean be fou in article by Hull and Hull [2]; the art ible 1 has tx ke 
a S/R e4-Nih SF 24-S SF S 
0 1. 38620 { ) 0314 L™ DSi7 
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~- 1, 23471 1516 1.» 1. 1000 28 02449 
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1.1° K (the Brillouin entropy itself is, of course, 

ependent of 6, but the values of // are listed in 
ms of w/7/A, which appears in the analytical 
yressions for the exact entropy We thus have 
0.135 and a,=—0.1227273. The derivatives are 
nin columns 3 and 5 


Now, in any experiment, 7° is determined by 
erimental conditions, and 6 is usually known 
rly accurately (It is not illogical to derive an 


irate value of 6 from entropy data, the calcula- 
n of which involves 4; it is essentially the method 
successive approximation These together fix 
value of a, and a tabulation of AS’? is readily 
tained by making use of table 5 and eq (6 
One then plots a curve of AS/R versus (S/R), 
sing eq (4)) to obtain the correction curve appro- 


BLE 5 Brillouin entropy, entropy correction, and first 
atives with re spect to af 4/7) for the standard value” 
af 0.1227 273 0.135/1.1), for varying values of B ull 
A 5/2 (see tert 
B=ull sd Sa/f Spif LSJ ASI 
s6020 ) ni74 0. 1200 
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priate to one’s own experimental conditions. 
any one demagnetization, S,/R is readily obtained 
from the /7/T value and table 1, and AS’P is read 
off from the correction curve. 

To illustrate the efficiency of the method, table 6 
lists values of AS/R calculated rigorously and by 


0.113043 


cause a 


the above approximation for the case a 
(6=0.26 deg, T=1.15° K) and for the 
0.133333 (6=0.28 deg, T=1.05° K Inasmuch as 
one 1s only interested in an accuracy to the fourth 
decimal place, or two significant figures, it may be 
seen that the method is very useful. 





TABLE 6 Comparison of entropy-correction values calculated 
exactly and by the approrimate method for the cases (a 
a=0.1138043 bh 0.26 deq, T=1.15°K and h a 
0.133333 (6/k=0.28 deq, T=1.05°K 

" $3333 a 43 

B=yll/d 

LS Sk Sf SI 

‘ wy x va Approx 

0 0. OORS 0. 0ON 0. 0063 0. 0083. 

7 0059 wise Oo4 OAT 

ols ool (OO 139 O14 
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ise Equilibrium Relations 


and Zirconia-Titania 
L. W. Coughanour, R. S. 


The systems CaO-TiO, and ZrO,-TiO 
and by the observation of fusion characteris 
in the CaO-TiQO, system was confirmed. One 
svstem. This compound, ZrO ,-TiO:, has bee 


the following parameters at room temperature: a 


solution development was observed in the 


svstem 


I. Introduction 


\ study of phase relationships in the systems CaQ- 
TiO, and ZrO,-TiO, has been conducted as a part of 
a program of fundamental studies of ceramic dielec- 
tries. Both of these systems have been investigated 
by previous workers. The results of these investiga- 
tions were rather divergent, however, so it was con- 
dered essential that the systems be reinvestigated 

lhe number of compounds proposed for the CaQO- 
TiO, svstem hes varied from one [1]? to three [2] 
evious studies of the ZrO,-TiQ, system |8, 4] 
wlieated the absence of any compounds in the 
em and the presence of partial solid solutions. 
The possibility of the stabilization of the tetragonal 
form of ZrO, by TiO, has also been suggested [5]. 
Brown and Duwez [6] from X-ray and thermal-expan- 
sion measurements proposed a phase diagram indi- 
cating the compound ZrO,-TiO,. Their diagram, 
although not supported by fusion-point 


Ss 


iit 


data, is 
similar to the one determined in the present Investi- 
gation 

The determination of the melting points of the 
compounds, and of the solidus and liquidus tempera- 
tures at various points across the two systems, have 
supplied data from which a revised equilibrium 
diagram has been constructed for the systems CaO- 
TiO, and ZrOQ>-TiQ,. 


2. Sample Preparation and Test Methods 


The following materials were used as the source of 
components in the samples for this study: CaO 
Kimer Amend Co. reagent-grade calcium carbonate 
of 99.5-percent purity; ZrO,—-dense ZrO, of nominal 
99-percent purity, obtained from the Titanium Alloy 
Manufacturing Co.; TiO,—highly purified titania of 
purity over 99.9-percent, obtained from the Tita- 
nium Division of the National Lead Co. 

Calculations of weight composition were made to 

0.01 pereent, no correction being made for the 
percentage purity of the raw materials. 

s study was sponsored by the Office cf Ordnance Research, Department 


Army 
ures in brackets indicate literature references at the end of this paper 


were studied by 


From the data obtained, an equilibrium diagram is suggested for each svstem 
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Roth, and V. A. DeProsse 


means of solid-state reactions 


ties. The existence of two binary compounds 


» compound was found to exist in the ZrQ)-LiO 
‘n found to have orthorhombic svmmetry with 

4.806 A, b=5.032 A, « 5AAT A Solid- 
ZrQ,-TiO, system, but not in the CaOQ-TiO 


The components, in sufficient quantities to give a 
10-¢ sample, were weighed to the nearest milligram. 
They were then mixed with a binder of either a 5- 
percent soluble-starch solution or a saturated 
solution of paraffin in carbon tetrachloride, and 1- 
in.-diam disks were formed at a pressure of 5,000 


’ 
i 


’ 
i 


lb/in?. The pressed samples were fired for 4 hr at 
1,266° C on platinum foil in an air atmosphere, 
using an electrically heated furnace wound with 
80 Pt-20 Rh wire. The disks were cooled and then 


ground, repressed, and refired for 4 hr at 1,350° C, 
again under oxidizing conditions. X-ray examina- 
tion of the samples showed, after these two heat 
treatments, that the solid-phase reactions were 
essentially complete with the exception of the forma- 
tion of 83CaO-2Ti0,. The composition representing 
this compound was a mixture of CaO and CaO-TiO, 
after these two calcinations. However, upon heat- 
ing to 1,560° C, or higher, the reaction to form the 
compound 3CaQ-2TiO, took place quite readily. 
Following this preliminary heat treatment the 
disks were ground, remixed with starch-solution 
binder, and disks about \ in. high were formed at 
15.000 |Ib/in’ in a \-in.-diam mold. From these 
disks, the specimens for solidus and liquidus determi- 
nations were ground in the form of small four-sided 


’ 
‘ 


pyramids, grooved on each side [7]. The tests 
were conducted in a thoria resistor-type furnace 
iS, 9]. The specimens were placed in the furnace on 


a disk of a platinum-rhodium alloy for temperatures 
below 1,900° C. For temperatures higher than this 
value, the samples were placed on an iridium disk. 
The metal disk was on a support of ThO,. The 
heating rate during the was about 4 deg 
C/min until the solidus temperature was approached, 
when it was reduced to about 2 deg C/min, or less. 
There was an oxidizing atmosphere in the furnace at 
all times. The temperature and fusion character- 
istics of the specimens were observed with an optical 
pyrometer calibrated at the National Bureau of 
Standards. In the determination of the solidus and 
liquidus temperatures, the first sign of liquid forma- 
tion, rounding of the corners of the test pyramid, 
the solidus temperature. The 


tests 


was recorded as 


37 








temperature of complete melting, as arrived at by 
observations through the optical pyrometer and by 
visual inspection of the cooled specimen, was re- 
corded as the liquidus temperature. In all cases, 
several had to be made to determine these 
two values 

The determination of such temperatures is subject 
to a number of sources or error. Among these are 
the slight deviations from blackbody conditions, the 
introduction of small quantities of impurities into 
the specimens in the forming and grinding opera- 
tions, the possibility of reduction of TiO, at high 
temperatures, and the inherent difficulty in’ the 
visual determination of the solidus and liquidus 
temperatures. In the present study, care was taken 
to keep each of these possible errors to a minimum. 
It is believed that the temperatures as recorded are 
accurate to + 10° C, 

After the fused specimens had cooled, they were 
analyzed either petrographically or by means of 
X-ray powder patterns, using a high-angle recording 
spectrometer, and Cu Ka radiation. All X-ray data 
are reported as observed, uncorrected for tempera- 
ture or adsorption characteristics. Similar analyses 
were performed on specimens that had been quenched 


‘ 


from temperatures up to 1,600° C. 


tests 


3. The CaO-TiO, System 


A considerable amount of work has been done on 
this system by previous investigators. The results, 
however, have been in rather serious disagreement, 
and the present study was undertaken to determine 
the correct equilibrium relations. 

Tanaka [1] decided that there is only one com- 
pound, CaO-TiO:, in the system and that CaO 
forms solid solutions with this compound. Hedvall 
and Anderson [10], by means of X-ray studies of 
samples fired to a maximum temperature of 1,135° C, 
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1. Suggested equilibrium diagram for the system 


CaO-TiOs. 


Ficure 








also concluded that CaO-TiQO, is the only compoun | 
in the system. H. von Wartenberg, et al. [2] pul - 
lished a liquidus curve for the system which indicate | 
the existence of three congruently melting com- 
pounds, 3CaQ-TiQ.,, 2CaO-TiO., and CaO-TiO 
Several other investigators [11-15] have indicated 
that the compounds existing in the system are 
3CaO0-2TiO, and CaO-TiO.. De Vries, et al. {15} 
proposed a diagram similar to the one here 
determined. 

The equilibrium diagram resulting from the 
present study is shown in figure 1. The fusion- 
behavior data from which the diagram is constructed 
are given in table 1. It is seen that there are two 
compounds in the system. These are CaO-TiQ, 
(perovskite), which melts congruently at about 
1,915° C, and 3CaO-.2TiO:, which melts incongru- 
ently at about 1,750°C. It should be noted that the 
compound 3CaQ-.2TiO, does not form readily at 
temperatures below 1,500° C. This may explain 
why it was not noted by some other investigators 
The X-ray powder diffraction pattern for the com- 
pound 3CaQ-2TiO, is given in table 2, but attempts 
to index this pattern have not, as vet, been success- 
ful. 

The absence of solid solutions in the diagram 
should be noted. There is considerable overlapping 
of some peaks in the diffraction patterns of the com- 
pounds CaO-TiO, and 3CaQO-2TiO,. This fact makes 
it rather difficult to unambiguously establish solid- 
solution relationships by X-ray methods. This 
difficulty is not present in the petrographic method 
of analysis, however, as their optical properties 
differ rather widely. Perovskite has very low 
birefringence and an index of refraction of about 
2.38. The optical properties of the compound 
3CaO-2TiO, may be tabulated as follows: Optical 


Fusion characteristics of compositions in the 
CaO-TiOz: system 


TABLe 1. 


Results of 


Results of petrographic 
heating c 


or X-ray examina- 
tion (specimens ex 
amined were not 
quenched and the 
phases observed may 
not be those in 
equilibrium at the 
liquidus) 


Composition 


Weight 


Mole 


Solidus —_ 


CaO TiO: CaO TiO: 


‘ r 
‘ ti 
4.76 05. 24 
10. 00 ©). 00 
2. 00 S80. 00 
25.00 75. 00 
33. 33 6. 67 


CaO. TiOe+ TiO, 
Do 
Do. 
Do 
Do 


Do 
Do 
CaO-TiO, 
Ca0-TiO.+ 
3CaO-27Ti0». 
Do. 


60. 00 
55.56 
5). 00 
47.00 


40. 00 
“4.4 
50. 00 
8.00 


55. 56 4.44 


Do. 
3CaO-2Ti0, 
3CaO-2TIO;+Cad, 

Do 

Do 

Do. 

Do. 

Do. 


43.00 
40.00 
36. 36 


57.00 
60.00 
6.4 
66. 67 
7.4% 
75. 00 
80. 00 


©. 00 10. 00 
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eter—biaxial | 2V ~80°; approximate in- 
of refraction—a=2.16, B=2.20, y=2.22 
values are in good agreement with those in 
rtial list of optical properties given by Fisk [14] 


2 Y-ray powder diffraction pattern for the compound 
3CaO.2TiO 
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4. The ZrO.-TiO. System 


\ liquidus curve for this system was published by 
von Wartenberg and Gurr [16]. They showed a 
eutectic point at about 80 mole percent TiO, and 
1,760° C. This is in exact agreement with the 
eutectic determined in the present work. Biissem, 
Schusterius, and Ungewiss [3] in an X-ray study of 
the system noted solid-solution development but no 
compound formation. Sowman and Andrews [4] 
published a phase diagram for the system Zro,-TiQ,. 
Their diagram was one of partial solid solution with 
no compounds present. They located the eutectic 
at 45 to 50 weight percent TiO, and 1,600° C. This 
diagram is not in agreement with the one determined 
in the present work. Later work by Sowman [5] 
on the ZrO,-TiO, system indicated to him that the 
|:1 composition, after heat treatment, consists of 
either a distorted form of tetragonal ZrO, or a new 
compound. Recently, Brown and Duwez [6] pro- 
posed a diagram indicating the existence of the 
compound ZrO,-TiO,. Their diagram is quite similar 


the one determined at the Bureau. 
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Figure 2 Suggested equilibrium diagram for the system 


ZrO.-TiO 


The equilibrium diagram for the ZrO,-TiQ, system, 
as determined in the present study, is shown in 
figure 2, and the fusion behavior-data from which the 
diagram was constructed are listed in table 3. It is 
seen that the compound ZrO,-TiO, melts incon- 
gruently at about 1,820° C. The eutectic pomt for 
the system is located at about 80 mole percent TiO, 
and 1,760°C. Rather extensive solid-solution devel- 
opment is noted at the higher temperatures. 

No attempt was made in this study to work out 
thoroughly the lower-temperature relationships in 
this system. Some observations on quenched speci- 
mens and slowly cooled specimens did, however, 
vield some information on this part of the system. 
X-ray analyses of specimens quenched from 1,500° C 
indicated that at this temperature about 12 percent 
of ZrO, will dissolve in TiQ,. This is in good agree- 
ment with the value given by Biissem, et al. (3). 


TABLE 3. Fusion characteristics of compositions in the 
/ 


ZrO 2-TiOs system 


Results of petrographic 
eating or X-ray examina 
tion (specimens ex 
amined were not 


Compositions 


Mok Weight quenched and the 
phases observed may 
Solidus Ua 1t be those in ec 
lidu ; not be th in equi 
Zrt rio ZrO rio mee librium at the liq 
an gi : . uidus 
4 ‘ ‘ 
oo. 00 10. 00 a3, 28 6. 72 Monoclinic ZrO 
75.00 25. &2.23 17. 77 1 aw Monoclinic ZrO 
ZrO, TiO 
6. 67 3.33 74. Sl 271.49 184 Ik 
ooo moO om wv. 19 1, S25 Do 
5 if 4.44 H5. 84 4.1 1, 820) Do 
f2.04 17.06 63. 44 1. 1, 825 1,870 Do 
0.00 50.00 #0. 71 19. 29 1, 820 1860) ZrO») TiO 
417.06 52.04 57. 82 42. 18 1, 700 L835) ZrOe TiOe+ TiO: ss 
44.44 » 55 23 “4.77 1, 775 1, 820 Do 
40.00 a0. 00 Ae 40.31 1. 70 1, 825 Do 
33.33 4. 67 3.54 Me. 46 1. 755 1.800 Do 
25. 00 75. 00 33. 05 66. 05 1, 765 1, 780 Do 
2. 00 80. 00 27. 83 72.17 1, 755 1, 760 Do 
16. 67 83.33 23. 57 7H. 48 1, 740 1, 780 Do 
10.00 90. 00 14.63 85. 37 1, 785 1, 805 TiO, ss 
5.00 95. 00 7.51 92. 49 1,810 1, 830 Do 











The mole composition 2ZrO,.:3TiO,, when quenched 
from 1,460° C, had only a small percentage of TiO, 
in it, and the parameters of the ZrO,-TiO, had 
changed, indicating some solid solution at this tem- 
perature. No solid solution of TiO, in ZrO,-TiO, 
could be detected in slowly cooled samples. 

Diffraction patterns of specimens high in ZrO, 
content indicated that TiO, is soluble in monoclinic 
ZrO, to the extent of 10 to 15 weight percent at room 
temperature. The ZrO, diffraction peaks in such 
samples had shifted to higher angles, indicating a 
smaller unit cell than that for pure ZrO,. This result 
is at variance with the findings of Brown and Duwez 
[6], who stated that the solubility of titania in mono- 
clinic ZrO, is probably very small. 


5. The Compound ZrO.-TiO. 


In view of the disagreement between the present 
authors and some previous workers concerning the 
existence of the compound ZrQO,-TiO,, some discus- 
sion of its properties is in order. First of all, it was 
noted that the 1:1 mixture formed a single phase 
when fired at 1,250° C for a sufficiently long period of 
time, or when fired at higher temperatures for shorter 
periods of time. There is a superficial resemblance 
between the X-ray pattern of this material and that 
of tetragonal ZrO,. However, the pattern of the 1:1 
fired mixture contained so many excess lines that it 
could not be indexed on the basis of tetragonal sym- 
metry. The patterns are compared in figure 3. It 
is evident that these patterns are for two entirely 
different materials. 

















Further, it was found that the 1:1 compositior 
when fused and cooled relatively slowly, gave crys 
tals large enough for petrographic examinatior 
Such examination showed the material to consist « 
a single crystalline phase having the following prop 
erties: Biaxial (—), 2V~80°; approximate indic« 
of refraction—a=2.33, B=2.38, y=2.41. 

These properties are not compatible with tetrag 
onal symmetry. It therefore must be assumed that 
the 1:1 mixture does not form stabilized tetragona! 
ZrO,, but does form a distinet compound. It is to bi 
noted that the X-ray pattern for the heated 1:1 mix 
ture is indexed on the basis of orthorhombic sym 
metry. The indexed pattern is listed in table 4 
The indexing of the powder pattern of ZrO,-TiO, was 
facilitated by the similarity of this pattern to that of 
tetragonal ZrO,. It was noted that the 002/200 
doublet of tetragonal ZrO, was split into three peaks 
indicating a lower symmetry for ZrO,-TiO,. When 
this triplet was assigned values of 002/020/200, the 
pattern was easily indexed and proved the compound 
ZrO,-TiO, to have orthorhombic symmetry. 

The theoretical possibility of the formation of the 
compound ZrO,-TiO, may be realized by noting the 
work of Dietzal and Tober [17]. Thev state that 
the tendency for compound formation between two 
oxides is dependent upon the relative difference in 
the field strengths of the two components involved 
The field strength may be represented as being 
proportional to the quantity za’, where the 
charge on the cation, and a ts the distance between 
the cation and anion. They noted that for most 
binary oxide systems A(z/a*) should be at least 0.30 
in order that compound formation may occur. 


Is 























40 50 6 
FURNACE COOLED ZrQ,: TiO, 











de 








cf lis Wer Lil ba, 
70 80 90 


202 








10 40 


20-—~ 


Fieure 3. 


sO 


TETRAGONAL ZrO, ss 


\-ray diffraction patterns for ZrO»-TiO, and tetrcgonal ZrO. 
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* R. I. is the intensity of diffraction peaks relative to the strongest peak 


Further, if compound formation does take place 
when A(2/a*) is less than 0.30, the compound must 
have an incongruent melting point. Furthermore, 
they indicate that A(z/a’) for the ZrO,-TiO, system 
is equal to 0.30. Thus, we see that the compound 
ZrO,TiO, is a possible occurrence according to this 
scheme, and that, if formed, it should have incongru- 





ent melting characteristics. This is exactly the 
situation observed in the present. study 

The compound when cooled slowly from a temper- 
ature of 1,200° to 1,600° C to room temperature has 
the following parameters: a=4.806A, 6=5.032A, 
e=5.447A. If the compound is quenched from 
1,200° to 1,500° C, the c-axis is found to have a 
larger value, whereas the a- and b-axes remain 
essentially the same. The parameters of the 
quenched compound are: a=4.802A, 6=5.034A, 
e=5.483A. If the compound is quenched from 
800° C, the parameters are the same as those in the 
slowly cooled specimen These results seem to 
indicate the existence of two polymorphic forms 
of ZrO,-TiO,, with the transition temperature 
between 800° and 1,200° C 

The only extinction rules which are observed for 
the low temperature form of ZrO,-TiO, are as 
follows; h0O, h=2n; OO, k=2n; 001, l=2n From 
these it must be concluded that this compound has 
a primitive lattice rather than a face-centered 
lattice, as is the case in the tetragonal, cubic, and 
monoclinic forms of ZrO, The high-temperature 
form of ZrO,-TiO, seems to have the additional 
extinction rule of: AOl, h+l=2n. This indicates 
that the high-temperature form of ZrO,-TiO, has 
a different symmetry than the low-temperature 
form. 


6. Summary 


The systems CaO-TiO, and ZrO,-TiO, were studied 
by observing solid-state reactions and fusion char- 
acteristics. The existence of two compounds in 
the CaQ-TiO, system was confirmed. They are 
CaO-TiO,, which melts congruently at about 1,915 
C, and 38CaQ-2TiO,, which melts incongruently at 
about 1,750° C. The two eutectic points were lo- 
cated as follows: 75 mole percent TiO, and 1,475° C, 
and about 23 mole percent TiO, and 1,725° C. No 
solid-solution formation was noted in this system 

One compound, ZrOQ,-TiO,, was found to occur in 
the ZrO,-TiO, system. This compound melts incon- 
gruently at about 1,820° C. The eutectic point for 
the system was located at about 80 mole percent 
TiO, and 1,760° C. Partial solid-solution develop- 
ment was noted in this system 

The X-ray powder diffraction pattern of the com- 
pound ZrQ,-TiQ, was indexed on the basis of ortho- 
rhombic symmetry. 
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Conrady’s Chromatic Condition ' 





Donald P. Feder 


Nearly fifty years ago Conrady published a simple and ingenious formula for determining 


the state of chromatic correction of an optical system 


matism”’, 


This method, called the ‘ D-d achro- 


depends upon the variation of the optical path as a function of wavelength. It 


can be used to evaluate both axial and lateral color, chromatic variation of spherical aberra- 


tion, and secondary chromatic aberration 
or as widely used as it deserves to be 


Unfortunately, 


this method is not as well known 


This paper proves the exact equation and discusses an extension of it which is approxi- 


mate. 


\ numerical example illustrates the accuracy obtainable with the approximation 


The principal applications are presented and formulas derived giving the relations between 
the “D-d method” and conventional means of expressing the chromatic aberration 


1. Introduction 


Conrady is among the best-known 
teachers of the art of optical design. His book 
\pplied opties and optical design’’ is familiar to 
most designers and is a bible to some. The “D-d 
method” for calculating chromatic aberration, how- 
ever, is not found in his book, and this may account 
for the fact that it is not as widely known as it 
rves to be. The method was published in two 
papers in the Monthly Notices of the Roval Astro- 
al Society nearly fifty years ago.” It was 
surprising to find recently that many good lens 
designers are unacquainted with the method.  Fur- 
thermore, even those who use it do not often 
apply it to off-axis imagery. 

lt will be shown in this paper, that after rays have 
been traced through the system to determine the 
monochromatic aberrations, the complete state of 
the chromatic correction can then be found with very 
little extra work and without the necessity of tracing 
rays at any other wavelength. In addition the D-d 
method gives the results of ray tracing at three wave- 
lengths. Because this method gives answers in the 
form of optical path differences, it is frequently 
easier to obtain a well corrected system by this 
method than by conventional methods. The use of 
the D-d method is recommended to designers who 
require a simple and elegant way for finding the 
chromatic aberration of an optical system over the 
entire field. 

The proof presented by Conrady is not entirely 
satisfactory, and later proofs, although more cumber- 
some, are less convincing. Some of them depend 
upon very complicated geometrical figures that are 
difficult to follow. 

The region of validity of the formula is not well 
defined. A recent paper, for example, stated that 
the chromatic variation of spherical aberration is not 
always given correctly by the formula. Conrady, 
himself, appeared to think that the secondary chro- 
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matic aberration * was exactly determined by the 
D-d method. It will be seen below, however, that 
any aberration determined by the first derivative 
with respect to wavelength is precisely given by this 
formula. Therefore, any chromatic aberration, in- 
cluding chromatic variation of spherical, is accurately 
given over a narrow wavelength region in the neigh- 
borhood of the median wavelength no matter what the 
monochromatic aberrations of the system may be For 
a somewhat wider wavelength region the method 
gives an approximate answer, which is ordinarily 
good enough to determine the secondary chromatic 
aberration of the system 


2. Derivation of Principal Equation 


Figure 1 shows a meridan section through a typical 
lens. The proof will be given for this case and 
immediately extended to any number of elements. 
A point source at O is emitting light that is refracted 
by the lens. For a particular wavelength , we 
consider a wave front in the image space. We may 
suppose that this wave front is located a considerable 
distance from the image of 0. We take a point, P, 
lying on the wave front, and draw the normal to the 
wave front. This normal coincides with a ray in the 
image space whose extension in the object space 
passes through 0. This ray is shown as a solid line 
in figure 1. In like manner for any wavelength A, 
there exists a unique ray passing from 0 to P. 
Another such ray is shown by the dotted line in 
figure 1. 

In order to find the chromatic aberration, we 
should like to know what happens to the wave front 
when the wavelength is varied. We first consider 
what happens at the point P?, and then, by letting 
P move over the wave front, obtain the total effect. 
We first mark off along the ray the distances 1), 
D,, Ds, Dy intercepted on the ray by the various 
media. The distance )), is the distance along the ray 
from the last surface to the point P. As the wave- 


? Because of strong objections by some people to the term secondary spectrum, 
the expression secondary chromatic aberration is used for this quantity throughout 
the paper 











Fieure 1 Veridian section of typical lens. 


Light from a point source at O traverse is producing an emergent 

e front at P in light of wavelength 1¢ shows such a ray with 
various media indicated by Dy, Dy The dotted 
from O to P for another wavelength A The wave front 


figure 


im the 


e distances 
melicates a ray going 


for A is not shown in the 
length changes, the distances 1)),, D))., Dy, Dy, change 
and so are to be considered functions of X. 

Let us define V(A) to be the optical-path length 
from O to P. Suppose, for convenience, that the 
object and image are in air; then 


VQ) =D, +N D+ ND Dy 


where .V, and NV, are the refractive indices of the 
corresponding media. Differentiating V with respect 
to \ gives the rate of change of the optical-path 
length from O to P. One obtains 


dV (dD, dD, | dD, dD) 
dy tdxr dx * dx . dr) 


, faN, 
i dr 


dNs 7 
D,+ an Dsx 
One sees that the derivative is divided naturally into 
two parts. In the first part (left-hand bracket), the 
refractive indices are held fixed. Therefore, this 
part results from the variation of the path of the ray 
and is zero by Fermat's principle. 

Fermat's principle states that the optical-path 
length between two fixed points possesses a stationary 
value with respect to arbitrary variations of the 
mechanical path. In this case we are considering a 
subset of all the possible variations in path from 0 
to P: namely, only such variations as are generated 
by varying \. Because of the manner in which this 
subset was chosen, it is clear that it contains, for 
light of any fixed wavelength \, the actual path 
traversed by the light. Furthermore, for light of 
this wavelength \, the actual path possesses a station- 
ary value. In mathematical language, 

dD, /dy+ NodD,/dyX+ NelD,/dy+ dD /dd=0. 
dV /dyk=D,dN,/ddX+ Dd N;/dX. 


It follows that 


| second term is more troublesome. 


This equation can be immediately generalized to any | 


number of elements giving 


dV dN 


dx a > dx D. (1) 


its*derivation restricts its application to axial obj; 
points or even to meridian rays. 


2.1. Significance of Principal Equation 


To understand the significance of eq (1) recall that 
for light of a fixed wavelength, the optical-path 
length, \’, from the object point to any point on the 
wave front is a constant. Suppose that dV/dX has 
been calculated for all points P on the wave front of 
wavelength \. For a sufficiently small increment in 
wavelength, Ad, it is approximately true that 
AV=(dV/d\)AX. A point P’ is now defined such that 
P’ lies on the ray through P, and such that the 
displacement of P’ from P is —AV. It is clear that 
the optical-path length V’ from O to P’ for light of 
wavelength \+ AA is equal to V. In this manner a 
set of points P’ is constructed, which constitutes a 
wave front for light of wavelength \+AX. By the 
nature of the construction, the displacement of this 
wave front from the original wave front is 

(dV d\)AX. This equation becomes exact only in 
the limit as AA. 

The above argument shows that if AV is a constant 
over the entire wave front, then the changed wave 
front is “parallel” to the original one, and there is no 
chromatic aberration in the image of 0. Stated more 
precisely, if dV dd, evaluated at A=, is the same at 
all points on the wave front, there is no chromatic 
aberration in the wavelength region around X». On 
the other hand, the variation of dV'/d\, as P is moved 
over the wave front, is a measure of the chromatic 
aberration. A more detailed study of this point is 
given later in this paper. 

The actual calculation of dV/dd from eq (1) is a 
simple matter. After the monochromatic errors are 
determined by tracing rays at a central wavelength, 
the values of the )’s are available as a byproduct 
The quantities dN/dd are available from the known 
properties of the glasses employed. It is a simple 
matter to form the indicated sum for several points 
on the wave front and to examine the chromatic 
aberration in the vicinity of the central wavelength 


2.2. The Approximate Equation 


Now one may ask how to determine the secondary 
chromatic aberration. For this it is necessary to 
calculate d?V/d\*. As eq (1) holds for all values of 
\, it may be differentiated, yielding 

CV /d¥Y = DEV /dN+S(dN/dd)(dD/dd). (2 
The first term of eq (2) is easy to obtain, but the 
Unfortunately, it 
is not zero, but computations on several widely 
different systems show that it must be considerably 
smaller than the term 2Dd* N/dd’. If this term, and 


all terms of higher order in AX, are neglected, then 
| one has 


This is the principal equation, and forms the basis | 


for the method. It should be noted that nothing in 
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dV 


dN (Ad)') 
4V= ar 


dV (Ad)? _ an 


nf dN 
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dx 
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AV=ZDAN. (3) 


tion (3) is an approximate equation, but it 
more information than eq (1) because it takes 
al account of all higher derivatives of N with 
et to ». Furthermore, it is much easier to 
than eq (1) because the glass catalogs list the 
es for discrete wavelengths, and in many cases 
he values of AN directly. On the other hand, 
ler to find dN/d x for use in eq (1) it is necessary 
one of the various dispersion formulas to the 
; data and to compute dN /d) for the wavelength 
; hich the rays were traced. Although this is a 
straightforward procedure, and has the advantage 
ing mathematically exact, the use of eq (3) not 
( entails less work but also gives an approximate 
of the secondary chromatic aberration 


3. Modification When Final Medium Is Not Air 


the case that the final image is not formed in 
then AVN’ rather than AV is the correct 
cpression for the displacement of the wave front 
Here N’ is the index of the final medium 
It must also be remembered that the sum SVAN 
should include the term D’AN’ for the last medium, 
where ))’ is the distance along the ray from the last 
refracting surface to the focal surface 


3. Application of Method 


\fter rays have been traced through the system to 
determine the monochromatic errors, it is necessary 


to determine the D’s in each medium for which 


AN#0. Fortunatelv, in the usual lens there are 
several ir spaces (W here AN=0), and these dis- 
appear from the summation in eq (3). In some 


ray-tracing schemes the values of the D’s are im- 
mediately available for each medium, and in 
others these quantities ere not directly given 
Even in the latter event, however, they can easily 
be computed from the ray-tracing data. A com- 
mon method for meridian rays employs the angle 
of incidence, 7, and the angle of inclination, U’, as 
coordinates of the ray. In this case one has 
D=r, cos I;—rz cos Ip+ (ty2—1r,4+-r2)cos Uy, 

where J; is the angle of refraction and r, the radius 
of the surface preceding the medium; /, is the angle 
of incidence, and r, is the radius of the surface fol- 
lowing the medium, and LU’), is the angle of inclination 
of the ray in the medium of thickness f,,». 

The values of the D's need not be obtained with 
high precision because each D is multiplied by a 
small number, which is rarely larger than 0.03. 
This is an important point from the standpoint of 
the practical computer. Five place accuracy in 
the D’s is more than sufficient for even the most 
precise work. 

\fter the values of the D’s have been computed by 
any suitable method, one forms for each traced 
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ray the sum SDAN=AV. 
for the region from F to C, then Np 
used. 


If one is achromatizing 
Ne for AN is 
Then AV is plotted against some variable 
that is a measure of the position of the ray in the 
aperture. For meridian rays this variable may be 
the coordinate of the ray in the entrance pupil. 
the tangent of the slope angle of the image ray, or 
any variable that locates the ray in the aperture 
For reasons of convenience, which will appear later, 
we shall use in this paper sin U’, where U” is the angle 
of inclination of the image ray 


3.1. Interpretation 
a. Axial Object Points 


If, for the axial bundle of rays, one plots the value 
of AV against sin U7, a curve that is symmetrical 
rbout the line sin (7=0 is obtained. The displace- 
ment of the curve is immaterial, and it is customary 
to subtract from each ordinate the value of AV at 
the origin. This value is SdAN, where d is the axiel 
distance between surfaces. One obtains 3(J)—d)AN, 
and this is plotted against sin 2’. This is the formu!a 
that gives rise to the popular name for the methed, 
the D—d method 

If the lens has no longitudinal chromatic aberra- 
tion, the curve will be a straight line (see fig. 2). 
If simple primary aberration is present, the curve 
will be approximately parabolic and will be concave 
upward in the case of overcorrection and concave 
downward for undercorrection. (In establishing a 
sign convention, AN is taken to be positive for all 
the media, and sin l’ is taken to be positive when 
the final ray slopes downward to the right. It 
should be noted that if AV is positive, the wave 
front for \+-AXd is behind the wave front for X.) A 
simple relation between the curvature of the parab- 
ola and the longitudinal aberration is derived in the 
final section. 

Usually the situation is complicated by the pres- 
ence of chromatic variation of spherical aberration. 
In many lenses initial undercorrection is followed 
by marginal overcorrection producing a zone. A 
case of this type is shown in figure 2. Good correc- 
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Figure 2. Arial fan. 


The ordinate for each curve is 2(D—d) AN As these curves are for axial points, 


they are necessarily symmetrical about the line sin U=0 
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Fieure 3. Oblique fan. 


The ordinate for each curve is Al Because the displacement is of no impor 
tance, each curve can be zeroed about the chief ray If the curve shows symmetry 
sbout the value of sin U corresponding to the chief ray, then only pure longi 
tudinal color is present A tilt indicates lateral color 


tion here makes ~(/—d)AN=0 at the margin. 
Then wave fronts that are tangent in the center of 
the aperture will cross at the margin. At a point 
usually near the 0.7 zone the wave fronts will be 
parallel and the corresponding rays coincident. 


b. Extra-Axial Object Points 


As mentioned above, for meridian rays one plots 
AV against sin (’, where L’ is the slope angle of the 
final ray in the image space and AV=SDAN. In 
this case, the symmetry present on axis is no longer 
inherent. Figure 3 depicts some of the character- 
istic types of chromatic aberration that may be 
present in oblique bundles. If the curve is a straight 
line inclined to the sin Ll” axis, then pure lateral 
aberration is present. One may think of the wave 
fronts as being inclined to one another. The lateral 
aberration is proportional to the slope of the line, 
and its numerical value can readily be obtained (see 
section 3.3). As before the displacement of the 
curve is immaterial, only the differences in AV being 
significant. If the curve is symmetrical about a line 
parallel to the AV axis through the chief ray, then 
pure longitudinal aberration is present. The longi- 
tudinal aberration for an off-axis object point may be 
quite different from that present on the axis. 

In the ordinary case the curve possesses neither 
symmetry nor straightness, and the aberration is a 
mixture of lateral and longitudinal chromatic 
aberration. 

In the final stages of an optical design, one may 
have curves of AV against sin L’ for several different 
object points. It is desirable to obtain the best cor- 


rection throughout the usable field. An advantage of 
this method is that the chromatic aberration is repre- 
sented in the same manner both on and off axis. 
The skew-ray aberration is also represented by AV, 
and in some cases it may be desirable to at least plot 













the sagittal ravs. The author believes, however, t! | 
the chromatic aberration of the skew rays will not je 
troublesome if the meridian rays are prope ly 
corrected. 


c. Secondary chromatic aberration 


To obtain an estimate of the secondary chroma ti 
aberration, one should recompute AV by using a 
different set of AN. Suppose, for example, that the 
central wavelength is sodium J and that rays have 
been traced in this color to determine the monochyo- 
matic aberrations. It is desired to correct the lens 
in the spectral region C to F. The procedure might 
be as follows: 

1. The sum 2(D-d)AN is made zero at the margin 
using Ne—Ne for AN and using the values of the )’s 
obtained from ray tracing in sodium light. 

2. Check 2(D—d)AN to make sure that it is satis- 
factory over the whole aperture. There is no point 
in worrying about the secondary chromatic aberra- 
tion if there is a large amount of chromatic variation 
of spherical aberration present. 

3. After =(D—d)AN is satisfactory over the 
aperture, it can be recomputed by using the value 
Ney—Np for AN and the same values for the D's as 
before. The value of this second sum is an indication 
of the secondary chromatic aberration. There will 
be some error owing to the fact that eq (3) involves 
an approximation, but this is not serious for most 
purposes. 

It should be recognized that the important use of the 
D-d method is to calculate the higher order chromatic 
aberrations. Before applying this method, a com- 
petent designer will already know the primary values 
of the longitudinal and lateral color and will have 
reduced these aberrations to reasonably small values 

If the secondary chromatic aberration is to be cor- 
rected, it must first be approximately corrected in the 
paraxial region. The simplest way to determine this 
is by tracing paraxial rays for the separate wave- 
lengths. After paraxial correction is made, it is 
necessary to apply the ))-d method as listed in steps 
1, 2, and 3 above to ascertain whether good correction 
has been achieved over the entire aperture. 


3.2. Tolerances 


One can set up a tolerance for the variation of AV 
on the basis of physical optics. Suppose that the 
image in J) light is sharp and that AV calculated from 
F to D does not vary more than one-fourth wave over 
the aperture. Then at the D focus the image in / 
light will be sensibly perfect within the meaning of 
the Rayleigh limit. Similarly, if the variation in 
AV calculated from PD to C is also less than one-fourth 
wave, then the image in C light observed at the )) 
focus will be sharp. Under such circumstances the 
chromatic aberration would be extremely well cor- 
rected over the entire region from F to C. On the 
other hand, the above tolerance is extremely strict 
and will seldom be met in practice. Usually one 
corrects the axial bundle of rays so that =(D—d)AN 
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»at the edge of the aperture, and this correction | angle of inclination that the final rav makes with 


inarily quite adequate. If there is a large zone | the optical axis.) Now one has 
ise of the presence of chromatic variation of 
cal aberration, then a basic change in the Ré R d(AV) 
, is necessary to reduce this aberration. If the AY ann t rw = and ds=Rdl. 
idarv chromatic aberration is to be corrected, 
sion for this correction must be made in the Then 
minary design. 
lenses covering an appreciable field, the curves R (d(AV)) | UAV 
several obliquities must be simultaneously cor- Ay=. 1) Ral ay W 
d, and this may lead to sacrificing best axial cos t Le 4 cos ( 
ction in order to gain at the edge of the field Finally 
D-d method is ideally suited to making such ti HAV 
sions because the chromatic aberration from Ay j ry 1) 
and extra-axial object points are both repre- d(sin 
sented by the same type of graph. 


The existence of eq (4) is the reason for plotting 
AV against sin (” rather then against some other 
aperture variable, such as tan [’ 

Equation (4) gives the displacement of the ray 
because of a change in A. One can plot AV against 


very many cases the tolerances are determined 
combination of experience and such overriding 
siderations as cost and complexity. 
































3.3. Relation Between Ray and Wave-Front sin (” for several traced rays. By measuring the 
Aberrations slope of the curve at any point, one finds the dis- 
placement between rays in the two wavelengths at 
For designers who are more familiar with ray that point. Using the formula, one can obtain the 
rcept methods for expressing chromatic aberra- | ray-intercept curve for any other wavelength after 
tion, it should be helpful to have a formula relating | the curve has been found in the central wavelength 
one to the other. Such a formula is derived below | In practice, however, it is not usually desirable to 
for meridian rays. do this because the single curve of AV against sin / 
One considers a wave front in light of wavelength | is easier to interpret than two ray-intercept curves 
\ located a distance ? from the focal plane (see fig. 4). | for different wavelengths. 
When the wavelength is changed by AX, the changed 
wave front makes an angle 6 with the original wave 3.4. Numerical Example 
front at some arbitrary point ?. It is clear that 
d(\V)/ds=6, where s is the are length measured Figure 5,a, shows a curve obtained in an actual 
along the wave front. Also s=U' because the | lens design. It represents the monochromatic aber- 
wave front is nearly spherical (Here l’ is the 
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(a) FOR CENTRAL COLOR 
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‘ |FOCAL — 00015 | 
\ wave PLANE 420 .380 .340 .300 .260 .220 
FRONT SINE OF FINAL ANGLE 
NN. Figure 5. 
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(a) The upper curve shows the meridian ray aberration at 60° obliquity from 
Fieure 4 the axis in D light. The ordinate H’ is the height of the ray on the focal plane, 
7 4 : the abscissa is the tangent of the final angle b) The lower curves show At 
for the same lens at the same obliquity. The solid curve is for the wave length 
full curve shows a wave front in light of wavelength A located a distance region F to D; the dotted curve from C to D 
the focus. A ray normal to this wave front at P makes an angle U with | These curves show the ray intercept curves at 60° in F light (solid line) and 
tic axis. A ray in light of wavelength 4+, makes an angle 6 with the | C light (dotted line) predicted by equation (4) from the curves in figure 5. The 
It is indicated by a dotted line through P crossing the focal plane at a | crosses show the actual results of ray tracing in F light (solid curve) and the circles 
¢ Sy from the first ray. | the results of ray tracing in C light (dotted curve 
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ration for J) light for an object point 60° from the 
optic axis. This lens was a complex system con- 
taining aspheric surfaces and was in an early stage 
of design. Figure 5,b, shows the curve of AV from 
F to D and the curve of AV from (to 2). Normally 
one would plot only the difference (that is, AV from 
F to C), but in this case the ray-intercept curves for 
all three wavelengths were desired. From the slope 
of these curves, the values of Ay have been calculated 


by eq (4). These differences were then added to 
the curve in figure 5.a, to obtain the curves in 


figure 6 and these curves checked by ray tracing. 
The crosses and circles show the position of the 
traced rays. The agreement is good except on the 
left end, where the aberration is large. The dis- 
crepancy is due to the fact that the wavelength 
interval used is not sufficiently small. 


3.5. Paraxial Correspondence 


One might ask the question how the curve for 
>(D-d)AN against aperture is related to the paraxial 
chromatic aberration. It would be possible to find 
this correspondence by an investigation of the limit- 
ing behavior of 2(/)-d)AN as the aperture is reduced 
to zero. Such an investigation, however, could only 
produce a well known formula giving the value of 
the longitudinal aberration in terms of the paraxial 
ray tracing data. This formula is reproduced below 
for reference. 

Instead of deriving this formula from an optical 
path difference standpoint however, we shall attack 
the problem in a different manner, that is, by 
assuming that the longitudinal aberration is known 
and determining =(J)-d)AN by using this informa- 
tion. Suppose there is a longitudinal difference in 
focus 6/ between the image positions for two wave- 
lengths. In the neighborhood of the axis each wave 
front is a sphere centered about its focus. Consider 
a fixed point on the axis, well removed from the 
focus, and the two spherical wave fronts passing 
through this point (see fig. 7). We ask the question, 
how does the distance between these spheres behave 
as a function of aperture? 

If z and y are the coordinates of a point on the 








Figure 7 


Two wave fronts in light of different colors in the neighborhood of the 
The longitudinal abberration is &. The wave front separation is A 


sphere of radius r, then z=}y’r~' is valid for a 
meridian section in the neighborhood of the axis 
where r is the radius. The z-axis coincides with the 
optical axis, and s=0 at the vertex of the sphere 
Differentiating with respect to r, one has Az { 
y’r-*Ar as the approximate displacement of the sphere 
of radius r+ Ar from the sphere of radius r. But 
Ar AV/N’, where AV=ZDAN and N’ the 
final index (see eq (3), section 2). Then AV=}N\’ 
y’r~* dl, since 6d/=Ar. Also sin U’=yr~', and so 


Is 


AV=3N'd sin’ UU. 5) 
Equation (5) shows that the curve of AV against 
sin (7 is a cirele or parabola for small values of 
and that the curvature is N’é/. This then is the 
paraxial correspondence desired. Another way of 
expressing this result is to consider AV as a power 
series in sin ©”. For an axial object point this series 
has the form AV=ce)+c. sin? U’+e, sint UU 
The odd powers of sin LU’ are absent because of 
symmetry. The first term ¢) is of no importance 
because it only represents a displacement of the 
reference point, which indicates that the wave fronts 
for the two wavelengths do not pass through the 
same axial point. This term ¢» is similar to the 
quantity SdAN, which it is customary to subtract 
from AV. 

The coefficient ce. depends upon the paraxial 
difference in focus, and by eq (5) e=—3.N’dl. The 
coefficients ¢s, ¢, ete. are related to the chromatic 
difference in spherical aberration. It might be 
emphasized that the D-d method gives the sum of 
this series precisely for small wavelength variations 
and approximately for finite differences in wave- 
length. 

In practice, the designer determines 4 either 
exactly from a paraxial ray trace, or approximately 
by a formula given below. Then AV is determined 
from the D—d method for a number of points 
throughout the aperture. (Frequently the values at 
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margin and the 0.7 zone are sufficient.) Then a 
ve can be drawn for AV versus sin U’ by making 
if the value of 4/ to draw the curve near the origin 
curves shown in figure 2 were obtained in this 
ner. 

e approximate formula for 4/ is 


él=a(N’u? 


a= > yi (AN. —wAN)). 


neanings of the above symbols, which all refer 
e paraxial ray, are as follows: 

s the paraxial inclination angle in the final 
im, 

s the height of incidence 

V/N’ at the refractive surface 


AN and AN’ are the differences in indices for two 
wavelengths for the media preceding and following 
the refracting surface 

The subscript k refers to the kth surface and the 
sum runs over ell the refracting surfaces in the 
system, the expression (N’u*)~' refers to the final 
medium 

Combining eq (6) with eq (5), one finds that 
AV=}a if one takes sin U’=u. This means that if 
one traces a paraxial ray through the system and 
computes a by eq (6), then AV along this ray is 
simply 4a. 


WasHINGTON, July 15, 1953. 





